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We study the collective excitations of polarized single-component quasi-two-dimensional dipolar fermions 
in an isotropic harmonic trap by solving the collisional Boltzmann-Vlasov (CBV) equation. We study the 
response to both monopole and quadrupole perturbations of the trap potential and investigate the character of 
excitations in each case. Simple analytic formulas are provided based on the linearized scaling ansatz and 
accurate numerical results are obtained by satisfying the first eight moments of the CBV equation. Except for 
the lowest lying monopole mode that exhibits a negligible damping in all of the studied cases, the quadrupole 
and the higher order monopole modes undergo a transition from the collisionless regime to a highly dissipative 
crossover regime and finally to the hydrodynamic regime upon increasing the dipolar interaction strength. For 
strong vertical confinements (2D limit), we predict the existence of a temperature window within which the 
characteristics of the collective modes become temperature independent. This behavior, which is a unique 
feature of the universal near-threshold dipole-dipole scatterings, persists as long as the scattering energies remain 
in the near-threshold regime. The predictions of this work are expected to be in the reach of current experiments. 



I. INTRODUCTION 

Dipolar quantum gases have been the subject of much inter- 
est and significant experimental and theoretical investigations 
in the recent years. The long-range anisotropic dipole-dipole 
interactions gives rise to novel phenomena and applications in 
these systems (see Ref. [1] and the references therein). In par- 
ticular, dipolar Bose-Einstein condensates (BECs) with mag- 
netic dipole-dipole interactions have been exhaustively stud- 
ied both theoretically and experimentally Q . The most re- 
cent experimental achievement along this line is the BEC of 
168 Er with a large magnetic dipole moment of 7 fis [3|. From 
an experimental point of view, the dipolar effects are much 
easier to observe in dipolar BECs compared to dipolar Fermi 
gases where the large energy scale set by Pauli exclusion re- 
quires larger dipoles for the interaction effects to become ap- 
preciable. Since electric dipole-dipole interactions are typi- 
cally stronger than magnetic ones, much of the recent exper- 
imental efforts have been focused on the realization of ultra- 
cold heternucleus bi-alkali molecules which have large per- 
manent electric dipole moments. 

One of the most important experimental achievements in 
this direction is the realization of a nearly quantum degener- 
ate gas of KRb molecules at JILA j4|. The experiments with 
other bi-alkali fermionic polar molecules such as LiCs |5] |5) 
are also making significant progress. At the same time, real- 
ization of degenerate fermionic atoms with strong magnetic 
dipolar interaction, such as 161 Dy [7], as well as microwave 
dressed Rydberg atoms (8J are other avenues toward realiza- 
tion of strongly interacting ultracold dipolar fermionic gases. 

An important experimental probe for the many-body 
physics of ultracold gases is the measurement of collective 
oscillations of trapped gases in response to perturbations of 
the trap potential. These oscillations constitute the low-lying 
collective excitations of these systems. The measurement of 
the frequency and damping of these oscillations can be uti- 
lized to understand the properties of the ground state and to 
extract important information such as the role self-energy cor- 
rections, the equilibrium equation of state, collisional relax- 



ation rates and kinetic coefficients. Moreover, the possibil- 
ity of carrying out extremely precise measurements of these 
quantities allows us to put our theoretical understanding of the 
system to the test. For instance, by measuring the frequency 
of the radial breathing mode for a two-component Fermi gas 
near the BEC-BCS crossover with a 10~ 3 accuracy level, the 
Innsbruck group could clearly verify the Quantum Monte- 
Carlo result for the unitary gas and invalidate the predictions 
of the BCS theory (9). Another good example is the recent 
measurement of the universal quantum viscosity of the uni- 
tary gas [ 10 1 that confirmed the theoretical T 3 / 2 scaling and 
also provided evidence for a conjecture on the lower bound 
for the viscosity over entropy ratio obtained using string the- 
ory methods ifTTIl . At the moment, the collective oscillations 
of trapped BECs [12 1 and two-component atomic gases with 
s-wave interactions in three dimensions |fl3l are both under- 
stood fairly well. Recently, the experimental and theoretical 
study of the 2D Fermi gas in the strongly interacting regime 
has also shown a significant progress Ifl4l - [l8l . 

Generally speaking, the low-lying collective excitations of 
an interacting system may be either described as collisionless 
(CL), hydrodynamical (HD), or in the crossover between these 
two limits. The CL limit is achieved when the the gas is either 
rarefied, or the interactions are negligibly weak or a certain 
dynamical symmetry forbids collisions. In this case, no dis- 
sipation occurs and the collective modes are undamped. The 
HD limit, on the other hand, is achieved either when the gas 
is in a superfluid state or in case of strongly interacting Fermi 
liquids, when the collision rate is much faster than the fre- 
quency of the collective modes so that a local equilibrium can 
be maintained [19] • In either case, the dynamics can be de- 
scribed well using simple HD equations in this limit, which 
are simply statements of conservation of mass, momentum 
and energy Ifl9ll20l . The collective modes are again dissipa- 
tionless in this limit. A realistic system, however, lies some- 
where between these two ideal limits, i.e. either the collision 
rate is not fast enough to maintain the local equilibrium or in 
the case of superfluid systems, the non-condensed component 
leads to collisions. An important aspect of understanding a 
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many-body system is to determine where it lies within this 
spectrum, both qualitatively and quantitatively. 

In this paper, we address this question for polarized 
single-component quasi-two-dimensional (quasi-2D) dipolar 
fermionic gases (see Fig. [2]i which has been the subject of 
much interest recently. In this setting, dipole-dipole interac- 
tions have a repulsive character and can be utilized to produce 
a strongly correlated Fermi liquid. Moreover, this particular 
configuration is found to be necessary in experiments with 
a wide range of bi-alkali polar molecules as a mechanism to 
suppress ultracold chemical reactions E). 

Several authors have already investigated certain aspects 
of this problem. In particular, Lima et al. have studied 
the collective oscillations in traps with various degrees 
of anistropy by assuming the validity of a hydrodynamic 
description GT1I221 while Sogo et al. studied the collisionless 
limit l23l . More recently, Abad et al. have studied both 
regimes separately and gave a comparison of the predictions 
of each 11241 . However, none of the mentioned works have 
given a quantitatively reliable condition for the validity of 
their approaches beyond simple order of magnitude analyses. 
Moreover, the intermediate regime in which one expects to 
observe the interesting physics of dissipation and collisional 
damping is not addressed in any of these works. 

In this paper, we make no prior assumption about where 
the system lies in the CL to HD spectrum. Instead, we use 
the framework of quantum kinetic equations, in particular, the 
collisional Boltzmann-Vlasov limit, which in principle allows 
one to study the dynamics in both limits in a unified way, in- 
cluding the crossover regime. The CL and HD limits naturally 
emergence when the right conditions are met. We evaluate 
the linear response of gas to monopole-like and quadrupole- 
like perturbations in the trap potential (x 2 + y 2 and x 2 — y 2 
respectively) and study the oscillation frequency and damp- 
ing of the generated excitations. We restrict our analysis to 
situations where the scattering energies lie well within the 
near-threshold regime so that Born approximation is applica- 
ble (25 26 1 . The collision integrals are treated without resort- 
ing to the usual relaxation time approximation. 

We carry out the calculations in two stages. First, we ne- 
glect the self-energy corrections to quasiparticle dispersions 
and utilize the widely used linearized scaling ansatz approx- 
imation ll27l to obtain a simple semi-analytic picture. In the 
next stage, we add mean-field corrections to quasiparticle 
dispersions and also extend the scaling ansatz approximation 
by satisfying all moments of the kinetic equation up to the 
eight order in order to obtain accurate numerical results. We 
find that both of these improvements result in significant 
quantitative corrections. Also, inclusion of higher moments 
also allows us to go beyond the study of lowest lying modes 
and to look at higher order modes as well. 

Before delving into the formalism and details, we find it 
useful to briefly summarize the main results of this work, 
some of which are unique and novel features of dipolar 
fermions in 2D. Without self-energy corrections, the scal- 
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FIG. 1. (Color online) Qualitative comparison of the dynamical 
regimes of quadrupole collective modes of 2D s-wave and dipo- 
lar fermions in harmonic traps, (al) and (bl) show the tempera- 
ture dependence of collisional relaxation rates of s-wave and dipolar 
fermions respectively. (a2) and (b2) show the dynamical regimes 
of quadrupole collective modes as a function of interaction strength 
and temperature. The asymptotics of ^~ wavc [ s due to 1151 . A?f 
is the trap Fermi momentum, ei2 is the 2D scattering length and 
Tt, — h 2 /(mkBCL^)- See Eqs. 1 1 3 1 and 25 for the definitions of the 
parameters appearing in (bl) and (b2). Please refer to the main text 
for details. 



ing ansatz analysis makes the well-known prediction of un- 
damped monopole oscillations at a fixed frequency of 2ojq, in- 
dependent of the interaction strength and temperature 112811291 . 
Here, loq is the in-plane (transverse) trap frequency. Taking 
mean-field corrections into account, we find that while the 
lowest lying monopole mode assumes only a negligibly small 
damping, its oscillation frequency will significantly increase 
from 2ujq due to the repulsive dipole-dipole interactions (see 
Fig. |6]l. The higher order monopole modes, however, are sig- 
nificantly affected by collisions: they go through a strongly 
dissipative crossover regime upon increasing the interaction 
strengths and finally approach the dissipationless HD regime. 

The quadrupole modes, including the lowest lying one, ex- 
hibit the same CL to HD transition as well. The oscillation 
frequency of the lowest lying quadrupole mode approaches 
V2wo m the HD limit, which is the universal, system- 
independent, frequency of quadrupole "surface" mode ll30l 
(see Fig. [8j. The emergence of surface mode is a clear in- 
dicator of hydrodynamics. 

We find simple semi-analytic results for quadrupole oscilla- 
tions using the linearized scaling ansatz and by dropping self- 
energy corrections. In this approximate picture, the frequency 
and damping of the quadrupole oscillations are controlled by a 
single dimensionless parameter, the collisional relaxation rate 



v c (Ref. to Sec. V 2 1. Small and large values of v c correspond 
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to collisionless and hydrodynamic regimes respectively. In 
the collision dominated regime, the viscosity sum rule yields 
v c ~ (P/r]s}trap, where P and r] s denote the local pressure 
and shear viscosity respectively |15|. By (. . .)t ra p, we imply 
averaging over the trap. Also, we identify v c ~ ^qt^ 1 in 
the thermal regime (T ^> Tp = y/2Nhu>o), where r c is the 
singleparticle collision time. 

For small T/Tp, we obtain v c ~ T 2 which is due to Pauli 
blocking. For large T/Tp, the behavior of v c depends on 
the degree of quasi-two-dimensinality (quantified by 77, see 
Eq. 13 1. In the strictly 2D limit [i] = 0), we find that v c 
reaches a plateau for T > Tp. The existence of this plateau, 
which is a unique feature of 2D dipolar fermions, results from 
the balance between rarefaction of the gas at higher temper- 
atures on one hand, and the growth of the dipolar scattering 
cross section on the other hand. The high temperature cut-off 
for this plateau is T^ip = S/(ma^g), where = mD 2 /ti 2 
is the "dipolar length". Here, m and D denote the mass and 
the dipole moment of a single particle. For T > Td; p , we find 

Fig. [T] shows a qualitative comparison between the behav- 
ior of quadrupole oscillations in 2D two-component fermions 
interacting via a s-wave Feshbach resonance (simply, s-wave 
fermions) and 2D dipolar fermions. The top and bottom pan- 
els show the temperature dependence of v c and the resulting 
dynamical regimes for the collective excitations as a func- 
tion of interaction parameters and temperature. The discussed 
regimes of v c for 2D dipolar fermions can be seen in panel 
(bl). It is worthy of mention that the temperature window 
in which v c is appreciably large is "universal" for 2D s-wave 
fermions. For 2D dipolar fermions, however, this window is 



amenable to experimental tuning (Ref. to Sec. V 2 1. 

We looked into the effect of mean-field corrections to quasi- 
particle dispersions and found that their inclusion yields sig- 
nificant correction in the quantum degenerate regime. This is 
again in contrast to the case of s-wave fermions where self- 
energy corrections are often found to have a negligible effect 
on the frequency of collective modes PP . 

Finally, going beyond the scaling ansatz and satisfying 
higher order moments of the CBV equation, we found that the 
scaling ansatz overestimates the collision rates in agreement 
to the findings of Ref. [31 1 in the context of s-wave fermions. 
We also found that the corrections to the energy of low-lying 
modes become negligible after 4th moments. The higher 
order modes were also briefly studied and we found that their 
behavior is qualitatively similar to the other modes. Finally, 
we discussed the experimental outlook of this work and 
gave predictions for the experiments with KRb. We found 
that although the HD regime is not achieved in the current 
experiments, there is a significant collisional damping rate 
which can be easily measured. 

This paper is organized as follows. In Sec. [TTJ we describe 
the model in detail and define the response functions. We 
briefly review quantum kinetic equation and the approxima- 
tions leading to the CBV equation and their validity condi- 
tions in Sec. sec:kinetic. We discuss the equilibrium state of 




FIG. 2. (Color online) A schematic picture of quasi-two-dimensional 
dipolar fermions in a shallow and isotropic in-plane (transverse) trap. 
Application of a strong dc field aligns the dipoles along the verti- 
cal axis (z). The quasi-two-dimensional limit corresponds to the 
[li/tmw,)] 1 ' 2 is much smaller than all of 
[h/ (mcJo)] 1 '' 2 , interparticle separation and the thermal 



situation where a. 
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de Broglie wavelength At = h/(2irmkBT) 
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CBV equation is described in Sec.[lV]and the variational cal- 
culation of the response functions using the method of mo- 
ments is discussed. The linearized scaling ansatz analysis in 
given in Sec. [VJ followed by the its extension to higher order 
moments and inclusion of mean-field corrections in Sec. |VI] 
Finally, we discuss the results in Sec.|VII| and the experimen- 
tal outlook of this work in Sec. I Villi Most of the technical 
details and tedious calculations are left to the Appendices. 



II. THE FORMALISM 

A. The Hamiltonian 

The Hamiltonian for trapped dipolar fermions prepared in 
a single hyperfine state and put in a strong polarizing dc field 
(electric for polar molecules, magnetic for atoms with perma- 
nent magnetic dipoles) can be written as: 



V2 

2m 



U, L> I d 3 rib\r)[-— + U™ p (r) j ^(r) 
+ /d 3 rdVV^(r-r')^(r)^t( I /)^ r ')^ r ) i (l) 



where: 



CCp(r) = \ mio 2 z z 2 + \mu 2 (x 2 + y 2 ), (2) 



is the axially symmetric trap potential and: 

D 2 



V d 3 P p (r) 



(3) 



the trapped gas in Sec. Ill The linear response theory of the 



We set h = 1 throughout this paper unless it appears explic- 
itly. A schematic picture of the system is shown in Fig. [2] 
We have assumed that the electric dipoles are polarized along 
the z-axis. Here, ip^> (r) is the fermion annihilation (creation) 
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operator is 3D space. In the limit uj z ^> luq, ep, ksT (where 
tF and T denote the Fermi energy and the temperature), the 
particles will be confined to the lowest transverse subband 
and we can reduce the above Hamiltonian to an effective two- 
dimensional model: 



H 2 » = J d 2 r^(r) (-^ + U™ p (r) + ^ ) ^ (r) 



+ J d J r dV V d f p (r - r') $(r) ^T (r /) ^ (r / } v , o(r) . (4) 

Here, r = (x, y) denote the 2D transverse coordinates and 
ipQ (r) denotes the fermion annihilation (creation) in the low- 
est subband. We have neglected the constant zero point energy 
hio z /2 of the lowest subband. C/ t 2 ° p (r) = mio 2 (x 2 + y 2 )/2 is 
the transverse part of the original trap potential and V 2 P(r) is 
the effective dipole-dipole interaction in the lowest subband: 



V d 2 »(r) = J dzdz'|0oW| 2 |^o(^)| 2 V d T p (r^-z'), (5) 

where 4>q(z) = e _z2 /( 2a ^/(- v /7ra z )5 is the vertical wave- 
function of particles in the lowest subband and a z = 
(mwz) -1 / 2 is the transverse oscillator length. The above in- 
tegration can be done analytically and we find: 



v d ? P M 



2n2al 



-2 



r 



Aa 2 J 



(6) 



where K n {x) is the modified Bessel function of the second 
kind. In the momentum space, we get: 



v d 2 ; » 



2nD 2 



qa z e 



^/ 2 Erfc 



qg z 

V2 



(7) 



We remark that in the long wavelength regime (qa z <C 1), 
the dominant contribution to VjP results from the repulsive 
side-by-side part of dipole-dipole interactions. We denote 
V$P = V, Vjg = V and U%P = U in the remainder of the 
paper for brevity. 

It is worthwhile to study the behavior of the effective 2D 
interaction in various limits. For qa z <C 1, one finds: 



V(?) 



AV2ttD 2 



3a, 



2irD 2 q + 0(q 2 ), 



whereas for qa z 3> 1, we get: 

2D 2 ^k A :! 



3a, 



1 



q 2 a 2 z 



+ 0(g- 4 a; 4 ) 



(8) 



(9) 



Apart from the constant term in Eq. (jH}, which is immaterial as 
long as we are concerned with a single hyperfine state, we no- 
tice an initial linear growth with respect to q which eventually 
saturates to a constant constant for q ~ l/a z . We shall see 



later that this linear growth results in an interesting behavior 
for the low-lying collective excitations. 

In real space, for small r/a z , one finds a behavior similar 
to the 2D Coulomb gas: 



V(r) 



D 2 



|-2- 7 -ln[r 2 /(8a 2 )] + 0(r 2 lnr)} 



(10) 

and for large r/a z , the r 3 dipole-dipole interaction is recov- 
ered: 



V(r)^D 2 /r 3 + 0(a z /r 4 ). 
It is useful to define a "dipolar length": 

mD 2 

ad 



h 2 ' 



(11) 



(12) 



which is the length scale associated to dipolar interactions, as 
well as the following dimensionless parameters: 



(13) 



where ao = [fi/(mwo)] 3 is the transverse oscillator length 
and TV is the number of trapped particles. A<j is a measure 
of dipolar interaction strength and is of the order of the typi- 
cal value of interaction energy over the kinetic energy in the 
quantum degenerate regime, r\ is a measure of "quasi-two- 
dimensionality" and is of the order of the vertical oscillator 
length a z divided by the zero temperature Thomas-Fermi ra- 
dius of the trapped gas. The strict 2D limit ui z —> oo corre- 
sponds to r\ = 0. 



B. Linear response theory 

A typical experiment for measuring the collective excita- 
tions of trapped particles is the following: the gas prepared 
in a thermal equilibrium state at t < 0~. At t > 0~, 
the system will be subject to a local perturbation, such as a 
kick or modulation of the trap potential. A certain observ- 
able will be then monitored either with an in situ or absorp- 
tion imaging technique. If the frequency and amplitude of 
the perturbing potential is weak compared to the correspond- 
ing microscopic scales, such an experiment can be theoreti- 
cally investigated within the linear response theory. Let us 
denote the perturbing potential and the observable as 8U(r, t) 
and O(r) respectively, and their corresponding second quan- 
tized operators are SU = J d 2 r ipl(r) 5U(r,t) ifj Q (r) and 
O = d 2 r V'o(r) O(r) ip (r). The usual linear response the- 
ory then yields: 



(0) t = J At' J d 2 rdV X l(ry ; t-0O(r)«/(r',0, 

(14) 
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where Xdd( r > r '> t—t') is the retarded density-density response 
function: 

X* (r,r';t-t') = -i6(t-t')Tr{p [n(r, t), n(r',t)]}, (15) 

where n(r,t) — ^l(r,t)tp Q (r,t) is the density operator and 
p is the initial density matrix. At this stage, one may choose 
a proper many-body approximation scheme and attempt to 
evaluate the response function using the diagram technique. 
However, the lack of translational symmetry due to the pres- 
ence of the trap potential makes this method complicated. In 
practice, one will have to make assumptions about separation 
of microscopic and macroscopic time and length scales in 
order to proceed. It is, however, much more transparent to 
acknowledge the existence of such a separation of scales from 
the outset and reduce the complicated evolution equations 
of the non-equilibrium Green's functions to quantum kinetic 
equations. One may then formulate and evaluate the linear 
response functions in the language of quantum kinetic equa- 
tions. We describe this method in the next section, where we 
also briefly review the quantum kinetic equations approach, 
and conclude this section by defining the response functions 
relevant to monopole and quadrupole oscillations. 

The monopole oscillations can be excited by choosing 
5U(r, t) = 8U m (r, t) = A(t) mw„r 2 , where A(t) is the tem- 
poral shape of the perturbation (e.g. a <5-function, a finite pulse 
or a periodic modulation). We choose A(t) = Aq ujQ 1 S(t~) 
for concreteness. Also, the linear response to any other pulse 
shape can be determined from the impulse response. Note 
that we have defined the monopole oscillations as the response 
of the trapped gas to a ~ r 2 perturbation. One may choose 
any other isotropic trap perturbation (such as r 4 , etc). How- 
ever, such choices are expected to excite higher order modes 
as well, and not necessarily the lowest lying ones. Here, the 
observable is the variation in the size of the cloud, f 2 — (f 2 )o- 
We define the "monopole response function" as: 



xAt) =Ao 1 niLj 9(t) ((f 2 



(r 2 )o) 



(16) 



Likewise, we define the quadrupole oscillations as the re- 
sponse of the trapped gas to 5U(r,t) = 5U q (r,t) = 
A(t) muj 2 (x 2 —y 2 ) and define the "quadrupole response func- 
tion" as: 

Xx 2- V 2 (t) = A^mujo 6{t) (x 2 - f) t . (17) 
Note that (x 2 — y 2 ) due to the isotropy of the trap. 



C. From quantum kinetic equations to the collsional 
BoItzmann-VIasov equation 



Quite generally, the dynamics of confined quantum 
gases can be formulated and studied using the formalism 
of non-equilibrium Green's functions, i.e. either by solv- 
ing Kadanoff-Baym equations using a physically relevant 
conserving approximation |[T9l or by using the Keldysh- 
Schwinger diagram technique. Such a formulation, however, 



is only necessary when the spatial and temporal scales of 
inhomogeneities (the trap and its perturbation) are compara- 
ble to the microscopic scales. In experiments dealing with 
a large number of particles N in a shallow trap, there is a 
natural separation of temporal and spatial scales between the 
microscopic (single particle) and macroscopic (collective) 
dynamics. Exploiting this fact, one can reduce the compli- 
cated Kadanoff-Baym equations to the somewhat simplified 
picture of quantum kinetic equations using the well-known 
procedure of gradient expansion lfl9l l32l . 

There exist several decent treatments of the subject of quan- 
tum kinetics in the literature and we refer the reader to the 
excellent pioneering monograph of Kadanoff and Baym lfT9l 
and Ref. l32l for details. However, for the purpose self- 
containedness and in order to clarify the invoked approxima- 
tions, we provide a very brief review of the basic elements of 
the kinetic theory. Our starting point is the general quantum 
kinetic equation for a system composed of a single species of 
fermions (i.e. a gas prepared in a single hyperfine state): 

[me(G- 1 ) + ,iG^]-[i^, ( ReG + } = G < S > -G > I] < , (18) 

where G+(p,w;r,i) = (uj - p 2 /(2m) - U(r,t) - 
and G^(p, cj; r, i) are the retarded and greater/lesser non- 
equilibrium Green's functions in the mixed Wigner coordi- 
nates [32]. U(r,t) denotes the external potential (i.e. the 
trap) and is assumed to vary on a scale much larger than the 
microscopic scales. E + (p, ui; r, t) and (p, u>; r, t) are the 
retarded and greater/lesser self-energies. In the mixed Wigner 
coordinates, (p, w) and (r, t) denote to the Fourier trans- 
formed microscopic coordinates and the slow macroscopic co- 
ordinates, respectively. [^4, B] denotes the generalized Pois- 
son's bracket defined as: 

[A,B]=d u) Ad t B-d t Ad u B 

- V P A- V r B + V r A- V P B. (19) 

It is generally understood that G + encodes the spectral prop- 
erties of the system (single particle states) while G < and G > 
contains the information about the statistics of particles and 
holes, respectively. Likewise, the real and imaginary parts of 
S + describe the renormalization of the singleparticle disper- 
sion and the spectrum broadening while and describe 
the collisional scattering -in and -out rates. In analogy to the 
equilibrium case, it is fruitful to introduce the local spectral 
function A(p,uj;r,t), Wigner's function /(p, cj;r, t) and 
spectral broadening T(p, to; r, t) (hereafter, we drop the com- 
mon arguments of the functions unless it is necessary), such 
that G< = iAf, A = i{G> - G<) = -2 3m(G+) 
and T = i(E> - £<) = -2 3m(S+). The ki- 
netic equations can be partially integrated to yield 



(<?+)" 



p 2 /(2m) - «Ke(X^ 



iT/2 [19|. This 



observation, along with one's choice of a many-body approx- 
imation that gives the self-energies as a functional of G < and 
G > , and finally the kinetic equation (Eq. 18 1 for either of G < 
or G > constitute a closed set of partial integro-differential 
equations for / and A whose solution describes the slow 
non-equilibrium dynamics of the system. For the case of 
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<I>-derivable many-body approximations, the kinetic equation 
obeys differential conservation laws for mass, momentum 
and energy currents. Such conservation laws are essential for 
formation and propagation of collective modes fl9l . 



leads to collisions) on the collisional side. The retarded self- 
energy in the HF approximation is instantaneous and is given 
by: 



Although the formalism of quantum kinetics is much sim- 
pler than a full non-equilibrium treatment, it is still extremely 
difficult to solve them in reality without resorting to fur- 
ther approximations. One useful approximation relevant for 
weakly interacting systems is the quasiparticle approxima- 
tion. The idea is that in the quantum degenerate regime, 
only the particle-hole excitations near the Fermi surface are 
responsible for the slow dynamics. The lifetime of such ex- 
citations, r _1 (pi?,ei?), is proportional to Tp/T 2 which can 
be very large. Thus, one may safely neglect the spectral 
broadening of the Green's functions appearing in the Pois- 
son brackets and approximate the spectral function as A w 
2tt5(lu — p 2 /(2m) — U — This approximation yields as 
ansatz for the greater/lesser Green's functions: 

G< p (p, w; r, T) = 2m Z p 6(uj- E p ) n(p; r, t), 

G> (p, w; r, T) = -2tti Z p 8(u - E p ) [1 - n(p; r, *)], 

(20) 

where E p is the (local) quasiparticle dispersion and is ob- 
tained by solving Lu — p 2 / (2m) — U (r, t) — S + (p, E p ; r, t) = 
0, and Z p = [1 - <9 w £+(p,w = E^r,*)]" 1 is the (lo- 
cal) quasiparticle residue. n(p;r, t) = /(p, E p ; r, t) is the 
quasiparticle occupation number. Plugging this ansatz into the 
kinetic equation, we obtain the collisional Boltzman-Vlasov 
(CBV) equation: 



£ + [rc](p;r,t)= / d 2 r 



, dV 



(2tt)' 



V(r - r') 



~<5 2 (r-r')V(p-p')Jn(p';r',i), (23) 

where V(r) and V(p) are the two-body interactions in the real 
and momentum space. Dealing with long-range interactions, 
we have included non-local contributions in the Hatree term. 
Such contributions are clearly beyond the first order gradient 
approximation but their inclusion may be necessary for suf- 
ficiently long-range interactions (it is exactly the presence of 
such non-local contributions in the Boltzmann-Vlasov equa- 
tion for the plasma that leads to plasmon modes and Landau 
damping). However, we will show momentarily that non-local 
contributions are negligible for the case of dipole-dipole in- 
teractions. Also, note that since E + has no w-dependence, the 
quasiparticle residue is 1. 

The collision integral in the Born approximation is given 
by EH: 



I An] = 



d 2 Pl dV dVi 



(2tt) 2 (2tt) 2 (2tt) 2 



(2iry5 2 {AP)(2ir)6(AE) 



x - \M\ 2 (1 - n)(l - n^n'n 1 ! - nn x (l - n')(l - n' x ) 



0_ 
Of 



m 



V r + V p S+[n] • V r - V r E+[n] ■ V p 



(24) 



where M = V(p — p') — V(p — pi) is the Born 
scattering amplitude, AP = P + Pi — p' — pi 



and AE 

EL 



Er, 



E„ 



E„ 



Note that 



p 2 /(2m) + U(r,t) + S+[n](p;r,i). We have 
V r U{r,t) ■ V p )n(p;r,i) = I c [n]. (21) also used the shorthand n = n(p;r,f), m = n(pi;r,f), etc. 



I c [n] is called the collision integral operator and is given by: 

I c [n] = -iZ p [(1 - n) S< + n E>] . (22) 

The CBV equation can be thought as a generalization of the 
usual Boltzmann transport equation of classical gases by (1) 
including Pauli exclusion effect in the collision integral, and 
(2) self-energy corrections of quasiparticle dispersions. A 
crucial observation made by Kadanoff and Baym is that the 
one may use different conserving many-body approximations 
for left hand (known as convective or dynamical) and the 
right hand (collisional) sides of the kinetic equation, without 
breaking the conservation laws. Intuitively, the dynamical 
and collisional contributions describe different physics and 
as long as each respect the conservation laws, the conserving 
property of the kinetic equation is preserved as a whole. 

The main goal of this work is to study the effect of inter- 
actions to the leading order in the interaction strength on both 
collisionless quasiparticle transport and elastic quasiparticle 
collisions. We use the self-consistent Hartree-Fock (HF) ap- 
proximation on the dynamical side and the Born approxima- 
tion (which is the lowest order ^-derivable approximation that 



We conclude this section by discussing the validity of the 
adopted approximations. Since we have described the in- 
teractions using the lowest order diagrams, the predictions 
are quantitatively reliable only as long as the system is in 
the weakly interacting regime, i.e. <C 1 (see Eq. 13 1. 
For dipolar interactions, this condition is equivalent to di- 
luteness Ji^D^d *C 1, where ri2D is the 2D density and 
is the dipolar length defined earlier (Eq. [T2| . Dealing with 
a Fermi liquid with essentially short-range interactions (i.e. 
J d 2 rV(r) < oo), the major many-body corrections such 
as the screening of interactions and in-medium T-matrix cor- 
rections are expected to change the predictions only quanti- 
tatively and the leave qualitative features intact even in the 
strongly interacting regime (Ad 3> 1). Therefore, while we 
acknowledge the limitations our approach, we allow ourself 
to extend our analysis to A<j ~ 0(1) as well. 

Aside from the many-body physics, the validity of Born 
approximation in describing two-body scatterings and neglect 
of multiple scatterings must also be assessed. The Born ap- 
proximation is valid when hv ^> Va, where v is the typical 
velocity of the scattering pairs in the center of mass frame 
and a is range of interactions. Identifying a with ad and 
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[m max(fcsT, ksTp)} 2 , this his condition implies: 



max(fc B T, k B T F ) < T d = 



(25) 



where we have defined the "dipolar temperature" Tdi p . 
This is precisely the condition for near-threshold scatterings. 
Ref. |25l has studied the 2D dipolar scatterings in detail, both 
in the near-threshold and semi-classical regimes. The study 
concludes that Born approximation predicts the correct scal- 
ing of the total scattering cross section with respect to the 
scattering energy provided that mvad/fr < 0.3. Inclusion of 
multiple scatterings, however, results in 0(1) quantitative cor- 
rections as one approaches the semiclassical regime. In this 
paper, we confine our analysis to near-threshold scatterings. 
Therefore, the quantitative validity of our results crucially re- 



lies on Eq. (25 1. Finally, we assume that the following scale 



separation holds: 



Tp < r, 



dip 



«0 

a d 



(26) 



so that we can allow ourselves to investigate both the quan- 
tum degenerate regime (T/Tp <C 1) and the thermal regime 
(T /Tp > 1) up to T ~ Tdip. We note that this condition 
is well satisfied in the current experiments with both polar 
molecules and atoms with permanent magnetic moments. 



IE. THE EQUILIBRIUM STATE 

The first step in the linear response analysis using the ki- 
netic equations is to determine the equilibrium distribution 
about which the perturbation analysis is carried out. Notice 
the analogy with the linear response analysis using the di- 
agram technique, where the first step is the evaluate of the 
equilibrium Green's functions. 

As mentioned earlier, we assume that the external potential 
U(y) = moJ^r 2 /2 is independent of time for t < 0~ and 
the system is assumed to have reached a thermal equilibrium 
state. It is easily to show that the CBV equation has a unique 
equilibrium solution given by: 



n (p; r) 



exp 



~ + E (p; r ) + \ muj 2 r 2 - (i 



where we have defined the useful shorthand dT — 
d 2 r d 2 p/(27r) 2 . In the case of harmonic traps, it is useful to 
define the following scales coordinates: 



Po' 



r ^[2A7(m Wo ) 2 ] 1/4 , 
Po = [2N(mu ) 2 ] 1/4 . 



(29) 



In the scaled coordinates, the equation for the particle number 
is / df n (p;f) = 1/2, where df = d 2 r d 2 p/(27r) 2 . The 
equilibrium distribution function also reads as: 



^o(p;r) 



exp 







■S (r-;p) - n 



(30) 



where f3 = V2Noj /{k B T), fl = fi/(y/2Noj ), and the di- 
mensionless and scaled self-energy functional is: 

S+[n](p;f,i) ^oj^ 1 ( df'[V2JVV[r (r-r')] 



- <5 2 (r - f')V[p (p - p')]J n(p'; f, t). (31) 

The motivation for the introduced dimensionless coordinates 
can be understood by investigating the non-interacting equi- 
librium solution at low temperatures. In this case, one can eas- 
ily find analytic solutions for the equilibrium density, Uq (f): 



n o 0) ( r ) = / d^pn (p;r) = log 



1 + e ^-f 2 /2) /(27T0). 

(32) 

Integrating over r, we obtain the following equation for the 
chemical potential: 

M 2 + y f 2 + 2 T 2 Li 2 [- exp(-A/f )] = 1, (33) 

where T — /3 _1 . At low temperatures, the above equation ad- 
mits the solution p, = 1 - tt 2 T' 2 /6 + C(e _ ^)//3 2 . The zero- 
temperature Thomas-Fermi radius of the cloud is easily ob- 
tained from Eq. (32 1, yielding fl^i? = [2%/2A7(tow )] 1/2 = 
\[2tq. Also, the Fermi momentum at the center of the trap 
[2%/2/Y(mw )] 1 / 2 = \/2po- We note that 



u (o) 
is given by p y F ' 



(27) 



where we have introduced the shorthand Eq = S + [no]. The 
above equation has to be solved self-consistenty along with 
the expression for the self-energy, Eq. (23 1. It is easily ver- 
ified that the above solution satisfies 7 c [no] = and at the 
same time, it solves the left hand side of the CBV equation. 
The global chemical potential fj, has to be found such that the 
equilibrium distribution function yields the correct number of 
trapped particles: 



N does not appear explicitly in the above expression, and the 
expressions look formally similar. Moreover, the equilibrium 
distribution function has almost a finite support of radius 0(1) 
in the scaled coordinates at low temperatures (beyond which 
it becomes exponentially small). 

Once we take the interactions into account, we can no 
longer obtain analytic solutions and will have to find the equi- 
librium distribution function numerically. It is useful to in- 
vestigate the effect of non-local Hartree energy (the first term 



in Eq. 37 1 before we move on. Carrying out the momentum 
integration, Hartree contribution of the self-energy can be ex- 
pressed just as a function of the density: 



drn (p;r) = N, 



(28) 



S+[n](r,0 = W( 7 1 



d 2 f'v / 2iVV(r f')n(f - f',t). (34) 
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Observing that the density function is only appreciably large 
in a region of size 0(1) in the scaled coordinates and the ap- 
pearance of r ~ N 1 / 4 in the argument of interaction poten- 
tial, the above integral is expected to only depend of the values 
of the density within a small region of size ~ TV" 1 / 4 about f. 
Assuming that n(r) is a smooth function, we may expand n 
to quadratic order about f , yielding: 



E+[n](r,t) « ujq 1 / d 2 f' \/2TV V(r f') \n(r,t) 



r' • Vn(f , t) + f' a f' p d a d n(f, t)/2 



(35) 



The first contribution is the usual local density approximation 
(LDA) expression: 



lLDA 



-'H L 



(f,t) = 



n{r,t) / d 2 f' V(r r') 



«n( f , t ). 



(36) 



The gradient term vanishes due to the isotropy of V(r). The 
quadratic term is dominated by the long-range behavior of 
V(r) assuming that the short-range part of V(r) is integrable 
(which is the case for dipolar interactions, see Eq. [T0| . Ob- 
serving that the Hessian matrix of the density is also 0(1) in 
the scaled coordinates, we easily find that the quadratic den- 
sity variations yield a correction that scales like N 1 / 2 ~ a / 4 for 
a potential with power-law tail V(r) ~ r~ a . For dipolar inter- 
actions, a = 3 and we find that the beyond LDA corrections 
scale like TV -1 / 4 and become irrelevant for large N. Note that 
if we were dealing with an electron gas (a — 1), such correc- 
tions would grow larger with N and had to be retained. This 
is the reason that one has to treat the Coulomb interactions 
in its full non-local from when studying the transport in plas- 
mas; on the same note, we remark that the physics of Landau 
damping is expected to be absent with dipolar fermions in the 
thermodynamic limit. In the remainder of this paper, we treat 
the Hartree potential in the LDA approximation and use the 
following local self-energy functional instead: 



e iWM(p;M) 



x n(p';f',t) 



v(o)-vMp-p') 



A„ 



ay 

(27T? 



«(|P - p'\,v)n(p';r',t). 



In the last last, we have defined: 



u(x, rf) = 2ttx Erfcx 



XT) 
75 



(37) 



(38) 



where Erfcx(a;) = e x Erfc(x). The dimensionless param- 
eters Ad and 77 were defined earlier (Eq. 13 1 Note that the 



dependence on N enters the equations only through these two 
parameters. 



We obtain the equilibrium distribution function using a sim- 
ple iterative numerical method as follows: at the initial step, 
we set S = and define the function n (ji) = n[E ,p J ] 
(i.e. the distribution function obtained using the self-energy 
So = and chemical potential jT). Keeping £0 fixed, we 
find no such that J dT 710(^0) = 1/2- To proceed from 
i'th step to (i + l)'th step, we set = E + [ni], define 

rij+i(j[t) = n[Sj+i,/2] and for a fixed £4+1, we find fti+i 
such that J dr' n,+i(//,+i) = 1/2. At the end of the step, we 
set rii + i — > (1 — A)rij + A n^+i, where < A < 1. The last 
step is to stabilize the iterative procedure and damp possible 
oscillations that prevent convergence. We found the above it- 
erative procedure to converge to the solution in less than ten 
steps within a relative error tolerance of 10 -8 . 




FIG. 3. (Color online) Equilibrium quasiparticle distribution func- 
tion of quasi-2D dipolar fermions for various temperatures and in- 
teractions strengths (uj z = 2tt x 23 kHz, ujq — 2tv x 36 Hz, 
TV = 2200). (a) T/T F = 0.1, \ d = 0, (b) T/T F = 0.1, X d = 1, 
(c) T/Tf = 0.5, \ d = 0, (d) T/T F = 0.5, A d = 1. Red and blue 
regions correspond to occupied and empty states respectively. 

Fig.[3]shows the equilibrium quasiparticle distribution func- 
tion as a function of p and f for several temperatures and inter- 
action strengths. As one expects, the presence of interactions, 
which are effectively repulsive, results in the expansion of the 
gas in the trap (compare panels a and b) and thermal fluctua- 
tions smear the Fermi surface (compare panels a and c). 

The equilibrium density is shown in Fig. [4^. The expo- 
nentially decaying tail of the density at higher temperatures 
and reduction of the density at the center of the trap at low 
temperatures due to repulsive interactions can be clearly 
seen. We also compare the LDA and non-local Hartree 
self-energy functions in Fig. [4j3 for various number of 
particles in the trap. The relative correction to the density is 
in the order of 10~ 3 for realistic number of trapped particles 
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FIG. 4. (Color online) Equilibrium quasiparticle density of quasi- 
2D dipolar fermions (oj z = 2ir x 23 kHz, ljo = 2-7T x 36 Hz), (a) 
dashed and solid lines correspond to the non-interacting (Ad = 0) 
and interacting (Ad = 1), blue (top) and red (bottom) lines corre- 
spond to T/T F = 0.1 and 0.5 respectively. In all cases, N = 2200. 
(b) A comparison between the LDA (solid lines) and non-local 
(dashed lines) Hartree self-energy functionals. From bottom to top, 
= 500, 1000, 2200, and 5000. A d = 1 and T/T F = 0.1 in all 
cases. The non-local corrections are clearly negligible and become 
smaller for larger choices of N. 



and as argued earlier, becomes smaller for larger system sizes. 

Knowing the equilibrium state, we can move on to the in- 
vestigation of the low-lying collective excitations about the 
equilibrium state. To this end, we discuss the linear response 
theory of the CBV equation in the next section. 



IV. ANALYSIS OF THE COLLECTIVE MODES: 
LINEAR RESPONSE THEORY OF THE COLLISIONAL 
BOLTZMANN-VLASOV EQUATION 

The linear response can be evaluated using the CBV equa- 
tion by introducing a perturbation to the external potential, 
linearizing the resulting equation about deviations from the 
global equilibrium state, 5n(p;f,t) = n(p;r,t) — n (p;f) 
and solving the resulting linear integro-differential equation. 
The merits of this approach compared to the diagram tech- 
nique is the possibility of obtaining approximate solutions us- 
ing variational methods. 

Since we are mostly concerned with low temperatures here, 
it is fruitful to introduce the following ansatz for Sn: 



Sn(p;r,t) = 9(t) A (p; f) $(p; r, t), 



(39) 



where Ao = dno/dp, = /3no(l — no). We remark that the 
above ansatz is not restrictive at the moment and since for 
T > 0, Ao has a unbounded support and any arbitrary de- 
viation from the equilibrium state can be represented with a 
proper choice of The only exception is T = where Ao 
restricts the deviations to the local Fermi surface, which is in 
fact a favorable feature. Since the low-lying collective excita- 
tions essentially depend on the particle-hole excitations about 
the Fermi surface and observing that the pre-factor Aq peaks 



about the local Fermi surface, we expect the solution of the 
linearized CBV equation to be representable with a smooth 
choice of <i> [20). As we shall see, this feature allows us to 
construct decent approximate solutions by choosing a linear 
combination of smooth functions as a variational ansatz for 
<£. Plugging this ansatz into the CBV equation, expanding 
to first order in $ and taking a Fourier transform in time, we 
obtain the following linear integral equation for $(p; f , u): 



iwA $ + g[<5>] 



- w 



<{n ,5U(r f,(j)}, (40) 



where {(f), ip} = V r </> • V p -0 — V p • V r ?/> is the ordinary Pois- 
son bracket and ui = uj/uj a . describes the collisionless 

self-consistent mean-field dynamics of quasiparticles: 



= A {$, U } + {n , E[A $]} 
= Ao{* + E[A $],7M, 



(41) 



where T-L — (p 2 + r 2 )/2 + £ - To get the second line, we 
have used the identity {no, A} = —Ao{7io, A} which can be 
easily proved by direct calculation and is valid for arbitrary A. 
The first term describes the dynamics in the equilibrium mean- 
field. The second term describes the the mean-field generated 
by the deviations and is a consequence of our self-consistent 
treatment. describes the collisional dynamics and reads 



d 2 pi d 2 p' d 2 pi 



2 J (2tt) 2 (2tt) 2 
x (2n)6(AE) \M\ 2 S{$} n n ,i(l - 



(2^) 2 



(2tt)V(AP) 
-n' A (42) 



where AE = U (p, r) +K (pi, f) - U (p', r) -H (pi, ?), 

AP = p + p 1 -p , -p;,M = m(V[po(p-p')]-Vbo(p- 
pi)J), and S[$] = $(p;r,u>) + $(pi;r,w) - $(p';r,w) - 
$(p^;r, to). Note that we have included Hartree-Fock self- 
energy corrections in the collisions which is required to be 
consistent with the inclusion of mean-field effects in the col- 
lisionless dynamics. Specializing to the case of dipole-dipole 
interactions, we get: 

\M\ 2 = A 2 K|p - p'|,n) - u(\p - pi|,n)] 2 . (43) 
Formally, the solution of Eq. ( |40| can be written as: 

i {n ,SU(r r,uj)} 



$ = - (-zwA 



(2JV)3 



(44) 



and the linear response can be determined using Eq. ( 39 1: 

(0) t = Jar J ^ e -^A (p;f)$(p;f, W + )O(p;r). 

(45) 

The difficulty is in inverting the operator appearing in the 



parenthesis in Eq. (44 1. Decent approximate solutions can 
however be found using a variational technique known as the 
method of moments. To this end, we restrict the solution space 



of Eq. ( 40 1 to a subspace spanned by a set of basis functions 
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(also known as moment functions) {</> Q (p; r)} and expand $ 
and SU in this basis: 

$(p;r,w) = ^$ Q (w)0 Q (p;f), 

a 

(2N)-hu(r Q f, u>) = J2 Su <*(") MP; f). (46) 



Plugging this ansatz into Eq. (40 1 and evaluating the moments 
of the resulting equation with respect to each of the basis func- 
tions, i.e. multiplying the sides of the CBV equation by each 
of the basis functions and integrating over r and p, we find a 
closed set of linear equations for the coefficient set {'fa}: 



" t«{^ Q )Sa(w) + (MtaM) W a {u) + $ Q (w)] 
+ f^{£[A O 0a],#o})$a(w) ~ ^Ja*a(w) = 0, (47) 

where we have defined the "Ao-average" as: 

(A(p;r))) = JdT A (p;r).4(p;r). (48) 



Summation over repeated indices is implied in Eq. (47 1. 
The matrix elements of the collision integral, ,f a j3 
\ df 4> a <y[4>/3] can be put in the following symmetric form 
using the symmetry properties of the collision integral kernel: 



'a/3 



d 2 f 



d 2 p d 2 Pl d 2 p' d 2 pi 



(2tt) 2 (2tt) 2 (2tt) 2 (2tt) 2 
x (2tt)5(AE) (2tt) 2 (5 2 (AP) \M\ 2 S[0 a ] S[^] 

x 7iono,i(l-n{,)(l-n' Q>1 ). (49) 

The first term on the second line of Eq. ( |47] i can be put in 
a more useful form using the identity <f>p{E[Ao<f> a ], Hq} = 
{^E[Ao0 o ],«o} - S[A o Q ]{^,H o }- Taking the Ao- 
average of both sides on this identity, we find that the first 
term on the left hand side vanishes. To see this, note that 
({if>,Ho}) = / df AoW-,Ho} = J df {A ^, Ho} for arbi- 
trary ip. The last equality holds since {A ,T-lo} = 0. Since 
Ao — > exponentially fast for large r or p, the divergence 
theorem implies that the last integral vanishes as long as i\> is 
exponentially bounded. Here, ip = <^E[An0 a ] which is in 



fact exponentially bounded. Finally, Eq. (47 1 can be put in the 
following matrix form: 



Hn 



l o )*(w) = -H *U(w), 



(50) 



where: 



(M) Q(3 = (<l> a <i>p), 

(H ) Q /3 = ((^{0/3,^0})), 

(Z) Q/3 = (S[Ao^]{0 a ,^o})), 

0c)a/3 = 



(51) 



and 3>(cj) and 5U(w) are the vectors with entries $ a (w) 
and <5v Q (w) respectively. If the observable 0(p;r) is also 
expressible in terms of the basis functions, 0(p;f) = 



2^ a O a (^ Q (p; r), then the linear response can be expressed 
as: 

(0) u = J dTOpfa A o $ Q (w+)0 Q 

= O t M*(w+). (52) 



Eqs. ( 50 1-( 52 1 are similar to the analysis given in Ref. OT1 for 
the case of s-wave fermions. Here, however, we have an addi- 
tional matrix Z that accounts for the self-energy corrections. 

It is useful to express the "evolution matrix", which we de- 
fine to be E = M _1 (Ho Z l c ), in its diagonal basis: 



iVftV~ 



(53) 



where ft is a diagonal matrix. Note that in general, E is a 
not a Hermitian operator and may have complex eigenvalues. 
Moreover, it is a non-normal matrix and therefore, its eigen- 
vectors are not orthogonal. As a side note, the same non- 
normality feature of the linearized BV equation for plasmas 
lead to Landau damping [33]. However, as we argued earlier, 
dipole-dipole interactions are not long-ranged enough to give 
rise to such effects. 



Using diagonal form of the evolution matrix, Eq. ( 50 1 can 
be expressed as: 



#(w) 



-iV 



1 



ft 



1 H q SV(uj). (54) 



The real and imaginary parts of ft determine the oscillation 
frequency and damping of the corresponding eigenmodes. 
Clearly, not all of the eigenmodes are expected to get excited 
in response to a given perturbation. This becomes particularly 
important when one is dealing with a large variational basis 
set. In such cases, as we will see later, the evolution matrix 
will have poles which are very close to each other on the com- 
plex frequency plane and it is not a priori clear which one(s) 
and in what proportion contribute to the response of the sys- 
tem. Using the linear response formalism, however, we don't 
have to deal with this question separately. Using Eqs. (52i 
and (|54|, we get: 



(0)u 



E r a (uj) 

a. 



r a {u) = -i[v T Moyv- 



l H ^UH] c 



(55) 



i.e. the residues r a can be expressed in terms of the known 
matrices. Note that for a Dirac delta perturbation in time, 
JU(w) in constant and independent of lj and so is r a (u>). 

Our goal here is to evaluate the linear responses accurately 
within the approximations made so far. In practice, the reli- 
ability of the approximate linear response functions obatined 
using the method of moments depends on one's choice of 
the basis functions. This choice can be motivated by the 
symmetries of the perturbing potential and the equilibrium 
state. Here, the trap potential is assumed to be isotropic 
and it is easy to see that [— icuA + S> — J? , S z ] = , 



where S z 



s: 



(r) 



Si p) , and S 



(r) 



i{xd y - yd x ) and 
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= i(p x d Py — PydpJ) are the rotation operators in the 
coordinate and momentum space respectively. Therefore, if 
SU lies in a certain eigenspace of S z , so will the solution 
of the linearized equation $ and one may choose the basis 
functions within the same eigenspace. Another symmetry 
which is preserved by the CBV equation is the reflection 
symmetry. More explicitly, defining the reflection operator 
as R x 4>(p x ,p y ;x, y) = <fi(—p x ,P v ] -x, y), it is easy to show 
that the linearized evolution operator commutes with R x as 
well. We will utilize these observations to define appropriate 
(and extensible) basis sets for monopole and quadmpole 
oscillations. 

Before we attempt to present accurate solutions which in- 
evitably requires heavy numerical calculations, we find it use- 
ful to make simple analytical predictions using a small basis 
set as first step. We use the scaling ansatz approach to find 
such a basis set and neglect self-energy corrections to sim- 
plify the calculations. We extend the basis set and include 
self-energy corrections afterwards and discuss the nature and 
importance of corrections that arise. 



V. LINEARIZED SCALING ANSATZ ANALYSIS 

The scaling ansatz provides a simple and intuitive picture 
of the collective excitations of confined gases. This method 
has been applied to various system in both isotropic and 
anisotropic traps, including Bose gases below and above T c , 
s-wave and dipolar fermions in the collisionless and hydrody- 
namics regimes ll2TH24l |27*1 . Here, we apply the method to 
the collisional Boltzman equation which as we shall see, al- 
lows us to study both CL and HD limits as well as transition 
from one regime to the other. 

In this method, one assumes that the non-equilibrium quasi- 
particle distribution function can be approximately described 
as a scaled copy of the equilibrium distribution: 

n sc .{p;r,t) = \ n [c/)^ 1 (p l - b.fjb^fjbi], 

(56) 

where hi and <f>i (i = x, y) are time-dependent scalings of po- 
sitions and temperature. The pre-factor is to ensure conserva- 
tion of particle number. The equilibrium solution corresponds 
to the choice b r = L, = d>~ = 



Here, (A, v) and (A, v) correspond to the isotropic (monopole) 
and anisotropic (quadrupole) scalings. Comparing the last 



equation with Eq. ( 39 1, we can recognize the first and sec- 
ond set of terms in the brackets as <& m o n and $q U ad, i-e. the 
variational ansatz that the linearized scaling ansatz provides 
for monopole and quadrupole modes respectively. 



The first term in Eq. ( 58 i, which is a consequence of the 
normalization prefactor of the scaling ansatz, requires fur- 
ther discussion. First of all, we note that this term in only 
non-vanishing in the monopole case. Since quadrupole oscil- 
lations are purely anisotropic, they do not violate conserva- 
tion of mass in the linear regime and therefore, no normaliza- 
tion results. The monopole oscillations as described by $ mon , 
however, violate the conservation of mass and the ansatz must 
be fixed with a counter term. The scaling ansatz fixes this de- 
fect with a uniform scaling of the distribution, leading to the 



first term in Eq. (58 1. 

We argue that such an ansatz is not a particularly good 
choice once collisions are taken into account and must be 
avoided since it may lead to unphysical conclusions. It is gen- 
erally understood that the non-equilbrium dynamics of degen- 
erate Fermi gases are governed by excitations near the Fermi 
surface while the fermions deep inside the Fermi sea remain in 
place due to their large excitation energy gap. A global rescal- 
ing of the quasiparticle distribution, i.e. a uniform rescaling of 
quasiparticle occupations irrespective of their energy gap im- 
plies mobilization of all fermions with the same probably, in- 
cluding those which are deep inside the Fermi sea (the density 
of states is constant in two dimensions). This is clearly an un- 
physical assumption and may lead to unphysically large col- 
lision rates. We note that it is well-known that the kinetic de- 
scription of the monopole mode in the absence of self-energy 
corrections, as we shall also show momentarily, must result in 
undamped oscillations due to conservation laws. 

To fixed this defect, we remove the global normalization 
factor and address the issue of mass conservation by simply 
allowing the chemical potential to vary instead. This amounts 
to adding a term ~ 6p,(t) dn /djl = A Sfi(t) to the ansatz, 
i.e. adding (j> = 1 to the monopole basis set. To summarize, 
we obtain: 



$mo„ = $Kt) + ci (*) r • p + c 2 (t) r 2 + c 3 (t) p 1 , (59) 



and: 



following re-parametrization of the scaling variables: 

b x (t) = 1 + X(t) + A(t), by(t) = 1 + X(t) - X(t), 
(j) x {t) =l + P(t) + v(t), cf>y(t) = l + 9(t)-p(t), (57) 



and expanding Eq. ( 56 1 to first order in A, A, v and v, we get 



Sn sc « -2(A + D)n + A [A r • p + Dp 2 + A f 2 ] 
+ A [A (xp x - yp y ) + v(p 2 x - p 2 y ) + A (x 2 - y 2 )] , (58) 

where in the last equation, we have neglected self-energy cor- 
rections for simplicity and used the non-interacting equilib- 
rium solution. Also, Aq = dno/djl = f3no(l — uq) as before. 



= 1. Introducing the $ quad = di(t) (xp x -yp y )+d 2 (t)(x 2 -y 2 )+d i {t){p 2 x -p 2 l ), 

(60) 

where <5/i(i), ct(t) and di(t) are to be determined. 

The determination of these unknown functions is usually 
done by plugging the ansatz into the kinetic equation, multi- 
plying the resulting equation by each of the basis function and 
integrating over the phase space to obtain a close set of differ- 
ential equations. This is equivalent to the formalism described 



in Sec. IV and we prefer to do it in the notation of this paper 



as a warm-up for the later sections where we extend the ba- 
sis set and include self-energy corrections. We finally note 



that the various terms appearing in Eqs. (59i and (60i can 
be simply understood in physical terms. In particular, r • p 
and xp x — yp y in $ mon and 'I'quad correspond to isotropic 
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and anisotropic macroscopic velocity fields, v n 



oc r and 



1. Monopole oscillations from the scaling ansatz 

Neglecting self-energy corrections, we have Y. = 0, Ho = 
(f 2 + p 2 )/2 and we easily obtain the following simple forms 
for M and H : 



Mi' 



1 







(((f-p) 2 )) 







\ 



and: 



iimon 

n 0,sc. 






\0 <p 4 



-2 -9 

r z p z 







(61) 



r 2 p 2 ) 



2(((f-p) 2 







(f-p) 






(62) 

where the basis is chosen in the same order as appears 



in Eq. (59i. The matrix elements of the collision integral 
vanish due to conservation of energy and momentum (see 
Eq.|49j and notice that S[l] = S[f 2 ] = 0, S[p 2 } = 2AE 
and S|r • p] = r ■ AP). Therefore, the oscillations will be 
undamped. 

It is possible find find analytic expressions for the A - 
averages appearing in the above matrices. However, using the 



relations 

//^2 -2 



and 



and «(r • p) 2 }) = 



^mon 



/o \ 
2 -2 
0-100 

\0 1 J 



(63) 



a result which is independent of temperature. The monopole 
excitation operator can be expressed is r 2 , which gives the 
"excitation vector" <5U = (0, 0, 1, 0) T in the monopole basis 



Using Eq. (50i, we 



(see the definition of <5U after Eq. 51 
finally find: 

$ mon (p; f , u) = [-2iw(f • p) + 2(f 2 - p 2 )} /{ui 2 - 4). 

(64) 

The frequency of oscillations is given by the poles of the de- 
nominator, ui nlon — ±2, which is a well-known result 11291 . 
We state without proof that extending the monopole basis has 
no effect on this result as long as self-energy corrections are 
neglected. In fact, it is a well-known fact that the full nonlin- 
ear Boltzmann equation (including collisions) admits an exact 
monopole solution with frequency 2cjo 1 29]. This is deeply 
related to the fact that the trap potential is harmonic and the 
particles have quadratic dispersion. Including interaction ef- 
fects or changing the trap potential will both result in violation 
of this result. 

We remark that besides the ui — ±2, the evolution matrix 
above admits two zero eigenvalues that correspond to eigen- 
vector $ ~ 1 and $ ~ f 2 + p 2 . Both of these eigenvec- 
tors correspond to unphysical excitations as they violate mass 
conservation. However, it is easy to see that both lie in the 
null space of H™°" at the same time. Therefore, in light of 
Eq. ( 54 1, they will never be excited regardless of one's choice 
of excitation operator 5XJ. 



2. Quadrupole oscillations from the scaling ansatz 

We find the following forms for M and Ho in the quadrupole 

basis: 



H^quad 



/ (f 2 p 2 )) 



p 2 } /2 (the first two of which are only valid for Harmonic and: 
traps), we find that they all factor out from the evolution ma- 
trix and we get: 



i_i quad 
n 0,sc. 



o 2«rV>; 

((f 4 )) o 
((P 4 )) o 



o \ 



P 4 ))J 







(65) 



(66) 



The order of basis functions is the same as it appars in 



Eq. (60 1. The only non-zero collision matrix element is J 2 ^, 
the rest of which vanish again due to conservation laws (see 
Eq. 



49 



and note that S[x 2 — y 2 } =0 and S[xp x 



VPy\ 



(xe x ~ ye y ) ■ AP). The collision integral can be expressed as 



follows using the results of Appendices C 4 and D 4 (in par- 
ticular, see Eq. |C20| >: 
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FIG. 5. (Color online) Frequency and damping of quadrupole oscillations of quasi-two-dimensional dipolar fermions in isotropic harmonic 
traps, (a) and (b): the frequency and damping of oscillations vs. v c respectively, (c) the damping rate of the overdamped component 
vs. v c . (d) the evolution of the damped oscillatory pole on the complex plane upon increasing v c in the range [0, 15]. (e) Q(T /Tp,rj) 

as a function of T/Tf for different values of r\ = (2iV)3 (cjo/^z) 3 - Q is related to the dimensionless collisional relaxation rate v c as 
v c — N(aa/ao) 2 Q(T/Tp,r]). The low temperature and high temperature asymptotes in the 2D limit are shown as blue and red dashed lines 
respectively. 



^3T ad = - 64^(2iV) ixjf 5 / p 5 dp 



271 



27T 11/ 



2?r Jo 2?r JO 



x sin 2 ((f) - 4>') 



d£ sin 7 £ cos £ / dv sin 5 v cos v 
Jo 

Xi Erfcx ( 2r]Xi \jTp) - X2 Erfcx ( 2rjX2 \l Tp 



1 



1 



cosh(p — fJ-/T) + cosh(p sin ^ sin 2v cos <p) cosh(p — /i/T) + cosh(/9 sin £sin2i/cos</>') 



(67) 



where xi = si n £ sin^|sin[(0 — 0')/2]| and X2 = 
sin^ sin;/ | cos[(0 — <j>')/2]\. The above integration can not 
be carried out analytically in general and requires a numeri- 
cal treatment. The analytical low T and high T asymptotics 
are given in Appendix [B| Note that the (dimensionless) non- 
interacting chemical potential p, is given implicitly by Eq. p3] l 
and depends only on the dimensionless temperature T, There- 
fore, except for the prefactor, the above integral is a universal 
function of T and rj. We define the "collisional relaxation 
rate" v c as: 



2^ _ {2N)^l 



Q(T,T)) 



= N[ ^ ) Q(T, V ). (68) 



The last equation also serves as the definition of the univer- 
sal function Q(T, rj). The quadrupole excitation operator is 
<5U quad = (0, 1, 0) T in this basis and finally, a simple calcu- 



lation similar to the monopole case yields: 

$ quad (p; r, uS) = [2u){v c - iui)(xp x - yp y ) + 2i(v c - iui) 
x (x 2 - f) + 2Co{pl - pi)] /D^ d {u, v c \ (69) 

where D quad (u), v c ), the "quadrupole characteristic equation" 
is: 



D quad (Q, v c ) = 0j{u) 2 - 4) + iv c (Cu 2 - 2). 



(70) 



The roots of D quad (o), v c yield the frequency and damping 
of the quadrupole mode. We note that Eq. (69 1, along 
with the characteristic equation given above, are "generic" 
results in the sense that one obtains the same expression 
for quadrupole oscillations independent of the specific form 
of interactions. For instance, Refs. Il28l and |34| obtain 
the same characteristic equation for s-wave fermions and 
a classical gas respectively. The model-specific details are 
encoded in the collisional relaxation rate v c . Therefore, it is 
worthwhile to study the generic features of the quadrupolar 



oscillations from Eq. ( 69 1 in terms of v c as a first step. We 
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return to the analysis of v c afterwards. 

Two important limits can be recognized for quadrupole os- 
cillations. The collisionless limit corresponds to v c — > and 
we find: 



lim $ quad (p;f, 



quad 
CL 1 



p;r,w) 



[ - 2iu{xpx ~ VPy) + 2(5 2 - f) ~ 2(PI-P 2 y )] /(u 2 - 4). 

(71) 

Notice the formal similarity to the monopole case. In this 
limit, we obtain undamped oscillations at ojQ^ d = 2ljq which 
correspond to the free motion of particles in the trap. In the 
limit of very fast collisions, v c — > oo, we find: 

lim $^(p;r,a;) = ^ u D ad (p;f,a;) 
= [ - 2iu(xpx ~ VPy) + 2(S 2 - f)} /{uj 2 - 2), (72) 

which describe undamped oscillations at a frequency 
w HD d = v^^o- This is the well-known quadrupolar 
"surface" mode which may also be obtained by solving hy- 
drodynamics equations for harmonically trapped gases l30l . 
The absence of damping despite that fact that the collision 
rate is very high can be understood by noticing that the 
mean free path of particles becomes much smaller than the 
system size in this limit and the only role of collisions is to 
maintain a local equilibrium state for each element of the 
gas. The entropy is locally extermal in this limit and there 
is no room for dissipation. This is essentially the physics 
of first sound. In contrast to the first sound, however, the 
surface modes only show up in confined gases. Although 
we have neglected self-energy corrections here, it is known 
that surface modes have universal frequencies since they 
correspond to divergence-less flows and are entirely driven by 
the trap restoring force 11301 . We will observe this universality 
in the later section, where we include self-energy corrections 
and obtain the same oscillation frequency in the HD limit. 

Except for the two ideal limits discussed so far, quadrupolar 
oscillations are generally damped for finite values of v c . This 
is due to the fact that the collisions are not fast enough to 
maintain the local equilibrium and thus lead to dissipation. 
The oscillation frequency and dampingrates can be found by 
analyzing the roots D quad (w, v c ). Fig.pa-c show the real and 
imaginary parts of the poles as a function of v c . In the limit 
v c <C 1, the poles are approximately at: 



± 2 



+ "T + l(D{v% (73) 



corresponding to a damped oscillatory mode at a frequency 
slightly less than 2ujq and a damping rate of ~ v c ujq/2. Ad- 
ditionally, there is an overdamped component with the same 
damping frequency to the leading order. In the other limit 
v c S> 1, we get: 

±(V2 + ^^j-~+0(u-% -iv c +iO(y- 1 ), (74) 



which describe a damped oscillatory mode at a frequency 
slightly higher than V^ujq and a damping rate of ~ v~ x ojq, 
accompanied by a (highly) overdamped component with 
a damping rate of ujqv c . Studying the residues of the 
overdamped poles, we find that the contribution of the this 
component is oc v 2 and oc v~ 2 to leading order in the CL and 
HD limits respectively. We associate the presence of such an 
overdamped component to the damping of initial excitations 
that lie far away from the local equilibrium. Finally, Fig. [5}i 
shows the evolution of the pole on the complex frequency 
plane upon increasing v c . It starts off at 2w i moves to the 
lower half plane and finally returns to the real axis at the 
hydrodynamic frequency \/2ui. 



We conclude this section by studying the behavior of 
Q(T, rj), which is the universal function that yields the 
collisional relaxation rate v c for dipole-dipole interactions 
(Eq. 68 1. In the collision dominated regime (i.e. v c ^> 1) 
where a viscous hydrodynamic description is admissible, the 
shear viscosity sum rule yiels v c as oJo(P/Vs)tr&p, where P, 
r] s and ujo are the local pressure, shear viscosity and the trap 
frequency respectively lfT51l . By (. . .}trap> we imply averag- 
ing over the trap. Also, in the classical regime (T 3> Tp), 
one finds v c ~ t^ 1 where r c is the typical time between two 
collisions 



We have calculated Q for several values of rj as a function 
of T by evaluating the integral given in Eq. ( 67 1 numerically. 
The results are shown in Fig. [5^. The asymptotic behavior 
of Q is investigated in Appendix [B] in the low and high tem- 
perature regimes in the 2D limit (?y = 0) and are shown on 
the same figure as red and blue dashed lines. We find that 
Q ~ T 2 for small T while it saturates to a constant value 
for large T. The low temperature T 2 scaling is related to 
the Pauli blocking effect, however, it is different from the 
case 2D s-wave fermions (and 2D electron gas), where one 
finds v c - T 2 \og(T/T F )~ 2 HUE]. This difference can 
be traced back to the fact that our system is spin polarized 
and the s-wave interaction channel is blocked. The logarith- 
mic enhancement of the shear viscosity (i.e. attenuation of 
v c ) originates from the logarithmic divergence of the s-wave 
scattering length in the near-threshold regime in 2D. We re- 
mark that the near-threshold cross section of all other scatter- 
ing channels remains bounded l26l . 

The high temperature plateau is a unique feature of near- 
threshold dipole-dipole scatterings in the 2D limit and its ex- 
istence can be understood in terms of the interplay between 
the temperature dependence of the scattering cross section and 
rarefaction of the gas. Provided that Tp <C T <C Tdi P , we can 
estimate the relaxation rate using a classical analysis by iden- 
tifying v c ~ 7 C = TjT 1 . The Born 2D scattering cross section 
can be estimated as <jb ~ 9 _1 |V((?)| 2 ~ qa 2 Erfcx 2 (ga z ), 
where q is the typical momentum of scattering particles and 
is ~ {inksT) 1 / 2 in the high temperature regime. The col- 
lision frequency is 7 C = t^ 1 ~ %^ m f p = hqna, where 
Z m fp = (ncr) -1 is the mean free path. The density at the center 
of the trap is jiq = muj 2 N/(2ttT) and decreases as ~ T _1 . 
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Combining these results, the collision rate amounts to: 



A. Variational basis set for monopole oscillations 



7c 



a 



Erfcx^ 



k B T 

hjj 7 



(T F < T < T, 



dip J 



(75) 

In the 2D limit, w 2 — > oo and we find j c — const (note that 
Erfcx(O) = 1). In other words, the growth of scattering cross 
section counteracts rarefaction of the gas to yield a constant 
collision rate. For finite vertical trap frequencies, the effective 
quasi-two-dimensional dipolar interaction weakens and we 
find that 7 C decays like ~ 1/T (note that Erfcx(x) ~ 1/x for 
large x). We remark that the single subband picture adopted 
here no longer holds true in the quasi-2D regime for large 
ksT / {huj z ) and one must include higher subbands into ac- 
count. We have shown in a previous work ||36l that all inter- 
subband interaction matrix elements have the same Erfcx fac- 
tor and therefore, we expect this qualitative behavior to remain 
unaffected. 

The plateau reached in the 2D limit relies crucially on the 
applicability of Born approximation. As mentioned earlier, 
the scatterings enter the semiclassical regime for T > Td ip 
(see Eq. 25 1 and Born approximation breaks down. In this 
regime, the total scattering cross section can be estimated 
using the Eikonal approximation l25l and one finds ersc ~ 
(ad/g) 1 / 2 . Repeating the same analysis with the semiclassi- 
cal cross section, we find: 



The generator of monopole oscillations, SU m ~ r 2 , belong 
to the zero angular momentum representation of S z . An arbi- 
trary function of such type can be expressed as f(p, r)[(x + 
iy)(Px — iPy)} n for arbitrary n G Z and f(p, r). Any smooth 
function of this type can be written as a power series expan- 
sion in r 2 , p 2 , r • p and £ = yp x — xp y . Observing that 
£ 2 = r 2 p 2 — (r • p) 2 , the most general basis for such func- 
tions can be constructed from the following two classes: 

4>t = 4>(m a ,n a ,k a ) = r 2m " p 2 ""(r-p) fc °, 
fc^K^W^^Air.pf". (77) 

Observing that R x <t>a = an d me f act mat me equi- 

librium state and the perturbations are reflection symmetric, 
we may discard {<fe}. We define as the "extended 

monopole basis" and drop the + superscript for brevity. To 
truncate the basis set, we keep all basis functions satisfying 
rn + n + k < M, where M is a positive integer which we call 
the order of the basis set. A first order basis set contains four 
elements, {1, r • p,p 2 , r 2 } and is equivalent to the linearized 
scaling ansatz discussed earlier. In general, a basis set of order 
M has (M+l)(M+2)(M+3)/6 elements. Expressions use- 
ful for numerical evaluation of the matrix elements of M, H , 
Y. and l c in the monopole basis can be found in Appendix [C] 



a 



knf 



(T > T dip ). (76) 



The qualitative behavior of v c for the full range of tempera- 
tures was shown earlier in Fig.[]Jl. 

So far, we have neglected self-energy corrections in the de- 
scription of the collective modes. We have also restricted our 
analysis to a variational calculation within a small basis set. 
In the next section, we extend our analysis to address both of 
these shortcomings. 



VI. EXTENDED BASIS ANALYSIS: 
THE EFFECT OF HIGHER ORDER MOMENTS AND 
SELF-ENERGY CORRECTIONS 



The general formalism described in Sec. IV allows one 
to account for self-energy corrections as well as obtaining a 
more accurate calculation of the response functions by ex- 
tending the variational basis set in a controlled way. Us- 
ing simple symmetry considerations, we introduce extensible 
polynomial-like variational basis sets relevant for describing 
monopole and quadrupole dynamics. Finite truncations of 
these basis sets allows one to satisfy all moments of the CBV 
equation up to the truncation order, which is an extention of 
our previous analysis. Since we are dealing with large basis 
sets and self-energy corrections at finite temperatures, resort- 
ing to numerical methods is inevitable and no simple analytic 
results are expected to be found. 



B. Variational basis set for quadrupole oscillations 

A quadrupolar function in two dimensions is a function 
that changes sign upon a simultaneous ir/2 rotation of both 
r and p about the z-axis. Such functions belong to the 
m z — ±2 representation of S z which can be expressed as 
f{p,r) e lM ^ r e lN ^ p , where M and N are two integers such 
that M — N = ±2, 4> r and 4> p are the angles r and p make 
with a fixed axis (we arbitrarily choose the a;-axis) and f(j>, r) 
is an arbitrary scalar function of p and r. One can identify 12 
classes of functions with such symmetry. Apart from an ar- 
bitrary scalar function f(p,r), the accompanying multipliers 
can be: 



— x 2 — y 2 , 



£2 — Px PyJ 



Vi = xy(ypx - xp v ), 



vt 



£J = xp x - yp y , 

PxPyiVPx -XPy), 



(VPX + XPy)(yP X - XPy), 



and: 



£1 = xy, £2 = PxPy, _ £ 3 = ypx + xp y , 

Vi = [yPx - xp y )(x 2 - y 2 ), % ee (ypx - xp y ){p 2 x -p 2 y ), 
% = (VP* - x Py)( x Px -yPy)- 

The functions with + and — superscript are even and odd 
eigenfunctions of the reflection operator R x , respectively. 
Like before, we can drop the odd class. Also, we find the 
following relations between these pre-factors: 

„2 C + 



2 ?? +=r 2 e 3 + -(r-p)? 1 + , 

2% + = (r-pK 2 + -p 2 £ 3 + , 
H =r 2 t+ -p 2 t+, 



(78) 



16 



using which we can drop the class of functions f{p,r)ru 
from the basis set. Since f(p, r) is assumed to be a smooth 
scalar function of p and r, in can be expanded in the monopole 
basis. Thus, in summary, we find that any smooth reflection 
symmetric quadrupolar function can be expanded in terms of 
\At §a } f° r i = 1) 2, 3 and a — (m, n, k), where m, n and 
k are non-negative integers and 0+ are the previously intro- 
duced monopole basis functions. We denote this basis set as 
the "extended quadrupole basis". We also remark that this 
basis set can still be reduced further in light of the relation 
2(r • p) g~ = P 2 £,t + r2 ^t' so tnat tne basis functions of the 
tv P e £2" r 2m p 2n (r • p) fc+1 can be written as a linear combina- 



tion of r im p 2n+2^ T . p )ft and £ 



+ „2m+2„2n 



p 2n (r ■ p) . Like 



before, we drop the + superscript for brevity in the remainder 
of the paper. An order- M truncation of the quadrupole basis 
set is the finite set that comprises all quadrupole basis func- 
tions satisfying k + m + n < M — 1. The first order basis set 
contains three elements, {x 2 — y 2 ,p 2 — p 2 ,xp x — yp y } and 
is equivalent to the linearized scaling ansatz discussed ear- 
lier. In general, a quadrupole basis set of order M contains 
M(M + 1)(2M + 7)/6 elements. Again, expressions useful 
for numerical calculation of the matrix elements of M, H , 51 
and l c in the quadrupole basis can be found in Appendix [D] 



C. Numerical results 

1. Preliminaries 

In this section, we discuss the numerical results obtained by 
evaluating the linear responses to monopole and quadrupole 
perturbations for various parameter. We vary and T/Tp 
in the range [0, 2] for fixed N = 2200. We study the 2D 
limit oj z — 00 as well as a quasi-2D case corresponding to 
the current experiments with KRb (ujq = 2ir x 36 Hz, lu z — 
2ir x 23 kHz |4|). This choice of parameters imply 77 ~ 0.322 
in the quasi-2D case. 

For each configuration, we performed the calculations 
within a 4th order basis set, comprising 35 and 50 basis func- 
tions for the monopole and quadrupole cases respectively, and 
satisfying all moments of the CBV equation up to 8th order. 
With the knowledge of the numerically obtained equilibrium 
solution (see Sec. [Ill), the matrix elements of M, Ho and Y. 
can be calculated with little computational effort using the ex- 
pressions provided in Appendices |C| and [P] 

The most computationally demanding part is the evalua- 
tion of the collision matrix elements. Although a consider- 
able number of them vanish due to either conservation laws 
or symmetries, a 4th order basis set yields 118 (monopole) 
and 307 (quadrupole) unique matrix elements each of which 
is a five-dimensional integral that has to be evaluated for each 
choice of Ad, r/ and T/Tp. This requires considerably more 
efforts and computation time compared to the simple scaling 
ansatz analysis we presented earlier, where only a single col- 
lision matrix element had to be dealt with. We calculated the 
collision matrix elements using the Monte-Carlo integration 
method with 5 x 10 s integration points. The estimate of the 
statistical error is ~ 10~ 3 (relative). 



We incorporate the self-energy corrections into the colli- 
sion integral within a local effective mass approximation (see 



Sec. C 4 1, which we found to be an excellent approximation 



in all cases. However, in order to assess the accuracy of this 
approximation and the consistency of the obtained results, we 
(1) we performed exact calculation of the collision integrals 
for a few representative parameter choices (without the effec- 
tive mass approximation), and (2) checked the satisfaction of 
conservation laws. We will discuss both of these consistency 
checks later. 

For the monopole case, we calculate the dimensionless 
spectral function A r 2 (w) defined as: 



A r 2(oj) = -(2N)-*3m\xAw)], 



(79) 



for a 5-kick in the potential (see Eq. [16). This quantity can 



be found using Eqs. (52 1 and ( |54] i by choosing the excitation 
and observation vectors as SU a (u) = O a = S ma , where m 
is the index that corresponds to the basis function <f> = r 2 . 
For the quadrupole case, we calculate the spectral function 
A x 2_ y 2(uj) defined as: 



A x 2_ y 2(u) = -(27V) i3m[x x 2- y 2{uj)] 



(80) 



(see Eq. [IT) . Likewise, this quantity can be evaluated by 
choosing the excitation and observation vectors as 8U a (u>) = 
O a = S qa , where q is the index that corresponds to the basis 



function 



£1 



y 2 . These spectral functions can be 



directly measured in the experiments in different ways (Ref. 
to Sec. [ynT) >. 

Although the evolution matrix has a large number of poles, 
some of which are clearly isolated and some may belong to 
branch lines, we found that only a few of them get excited and 
contribute to the response. Many of such poles lie inside the 
null space of Ho, are unphysical and do not get excited (see 
the discussion at the end of Sec. |V l| i. In all cases, we found 
that the spectral functions can be reproduced accurately by a 
fit function with two simple poles in the lower half plane: 
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A* 
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-Q-iT 



fl — iT uj — iT' 



(81) 

corresponding to damped oscillations and overdamped com- 
ponents. Such a fit function extracts the most important in- 
formation from the numerically obtained spectral functions. 
Moreover, this method allows us to present the results in clear 
and concise way. 



D. Monopole oscillations 



As mentioned earlier in Sec. |V 1| in the absence of 
self-energy corrections, the CBV equation for harmonically 
trapped gases admits an exact solution corresponding to a 
scaling velocity field v ~ r which has a fixed oscillation 
frequency of 2wo with no damping, independent of the in- 
teraction strength and temperature. This is due to fact that 
the Boltzmann equation admits a closed set of equations for 
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FIG. 6. (Color online) The oscillation frequency and the damping 
(inset) of the monopole mode extracted from the numerically ob- 
tained spectral functions using a 4th order basis set (including self- 
energy corrections). The colored and grayscale (upper and lower) 
graphs correspond to an ideal 2D system (77 = 0) and a quasi-2D 
system (77 ~ 0.322) respectively. Red and blue graphs correspond 
to high and low temperatures respectively. In all cases, N = 2200. 
The inset plot shows the damping rate in the 2D case (77 = 0). 




9te w/cdo 

FIG. 7. (Color online) (left) A typical picture of the poles of the 
evolution matrix (T/T F = 0.45, A d = 2, N = 2200 and 77 = 
0). (right) the mass currents associated to the indicated poles. The 
chosen poles correspond to the three modes with lowest energy that 
survive in the HD regime. 



the moments of r 2 , p 2 and r ■ p, all of which are immune to 
collisions due to conservation laws. Taking self-energy cor- 
rections into account, the quasiparticle dispersion relations no 
longer remain quadratic and one finds that this simple chain 
of moment equations can not be closed anymore. In particu- 
lar, contributions from higher order moments, many of which 
are strongly influenced by the collisions, become important. 
Therefore, we expect the monopole oscillations to be damped 
to a certain degree. 

Fig. [6] shows the frequency and damping of the monopole 
oscillations extracted from the numerically obtained spectral 
functions. The colored and grayscale (top and bottom) plots 
show correspond to the 2D limit (77 = 0) and the quasi-2D 
example (7/ ~ 0.322). The repulsive dipole-dipole interaction 
clearly results in a significant increase in the oscillation fre- 
quency. Also, as one expects, deviations from the 2D limit 
leads to a weaker repulsive effective interaction and thus, a 
smaller increase in the frequency of collective modes. 

We find that the dominant contribution to the response re- 
sults from a single isolated pole of the evolution matrix, which 
is the one that has the lowest energy. The relative residues of 
the other poles were found to be of the order of ~ 10~ 4 in all 
cases. 

The most interesting finding is that this mode exhibits a 
very small damping, 7, non < 10 _3 u;o> in all of the studied 
cases (see the inset plot of Fig. [6]) despite the presence of 
significant self-energy corrections. In Sec. VII we discuss 
the possibility that the smallness of damping could be a result 
of the effective mass approximation adopted in evaluating the 
collision integrals and confirm that even an exact treatment of 
self-energy corrections does not change this finding apprecia- 
bly. 



We remark that the mode which makes the dominant con- 
tribution to the linear response is the one that lies at the bot- 
tom of an infinite hierarchy of possible monopole oscillations. 
The reason that it is the only mode that gets excited is most 
likely a consequence of harmonic confinement. Inclusion of 
higher order moments not only yields a more accurate calcu- 
lation of the frequency of this mode, but also it allows one to 
investigate higher order modes. Fig.[7]shows a typical picture 
of the poles of the evolution matrix, along with plots of the 
mass current associated to three indicated low-lying modes. 
In contrast to the scaling mode (shown as a in the figure), the 
two other modes (b and c) have a significant damping rate. In 
the experiments, these modes can be excited by non-harmonic 
perturbations in the trap potential, such as ~ r A . 

In the absence of interactions, all of the poles lie on the real 
frequency axis at discrete locations 2nui , n € Z (here, up to 
\n\ = 4 due to the 4th order truncation of the basis set). Each 
discrete frequency is multiply degenerate (the collisionless 
Boltzmann equation, i.e. the Liouville equation, has infinitely 
many degenerate discrete poles in case of harmonic external 
potential). Upon increasing interactions, the poles spread to 
the lower half complex frequency plane, signaling the tran- 
sition to the dissipative CL-HD crossover regime. Upon fur- 
ther increment of the interactions, most of the poles diverge to 
— ioo while a few migrate back to the real axis and form the 
discrete hydrodynamic modes. The absence of damping in 



this limit, as discussed in Sec. V 2 is due to the emergence of 
local equilibrium. The three chosen poles in Fig.|7]correspond 
to isolated poles that survive in HD regime. 
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FIG. 8. (Color online) Evolution of the quadrupole oscillations from collisionless (CL) to hydrodynamic (HD) regime upon increasing the 
interaction strength (left to right). In all cases, T/Tf = 0.45 and r\ — (uj z = oo). The top row shows the spectral function and the bottom 
row shows the location of the poles of the evolution matrix on the complex plane. The pole shown as red is the pole that makes the dominant 
contribution to the response, (al) and (a2): Ad = 0.1, (bl) and (b2): Ad = 0.4, (cl) and (c2): Ad = 2. See Fig.|9]for a plot of the mass currents 
associated to the encircled poles. Refer to Sec.|VIIl|for a discussion on the experimental methods for measuring the spectral functions. 
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FIG. 9. (Color online) The mass current associated to the three 
modes marked in Fig. [8j;2. (a) is the lowest lying mode, known as 
the surface mode, characterized by the velocity field v ~ xe x — ye y . 
(b) and (c) are the next two modes. The nodal structure of the mass 
current is clearly noticeable. These modes constitute the three lowest 
lying HD modes. 



E. Quadrupole oscillations 

A typical scenario for quadrupole response is shown in 
Fig. [8] The top and bottom rows show the quadrupole spec- 
tral function and the location of the poles on the complex 
frequency plane respectively. The mode that has the largest 
residue is marked as red. For small interactions (\d <C 1, 
Fig. [8^1-2), the spectral function is sharply peak around 2ojq 
and the poles of the evolution matrix lie very close to the real 
axis. Upon increasing the interactions, the poles spread to the 
lower half complex frequency plane, indicating the entrance to 
the dissipative CL-HD crossover regime. The spectral func- 
tion is significantly broadened (see Fig. [HJ^l) in this regime. 
For stronger interactions, the local equilibrium picture starts 
to emerge, indicated by a reduction in damping. Fig. [8};2 
clearly shows a sharply peaked spectral function near \/2woin 
the strongly interacting regime. This is exactly the universal 
frequency of the hydrodynamic quadrupole surface mode dis- 



cussed earlier. 

Similar to the monopole case, we find that quadrupole per- 
turbations of the trap potential predominantly excite the low- 
est lying mode. Here, we find a small contribution from a 
few overdamped modes as well, specially in the crossover 
regime. This is in agreement to the result found from the scal- 
ing ansatz. Fig. [9] shows the mass current associated to the 
three lowest lying modes marked in Fig. [8};2. These modes 
are found to be the ones that survive in the strongly interact- 
ing regime and constitute the lowest lying HD modes. Such 
higher order modes may be excited by applying anharmonic 
perturbation to the trap potential, such as ~ r 2 (x 2 — y 2 ). 



Figs. 10 and 11 show the frequency and damping rate of 



the quadrupole oscillations obtained from the two-mode fit to 
the quadrupole spectral function. The result obtained from the 
scaling ansatz analysis presented earlier is also shown as thin 
black lines for reference. 

Clearly, the corrections are significant. In the low temper- 
ature regime, the self-energy corrections are dominant, yield- 
ing a a Ad correction to the frequencies (see the rightmost 
plot on the top panel of Fig. 
corrections are only oc 
(see Eq 



10 1. Note that the collisional 



73 



in the weakly interacting regime 
and note that v c oc A^). The corrections result- 
ing from the inclusion of higher order moments can also be 
seen in the same figure. Generally, scaling ansatz without 
self-energy corrections overestimates the collision rates and 
predicts the crossover to the hydrodynamic regime too early 
(notice that the peak of the damping rate occurs earlier for the 
scaling ansatz). At lower temperatures, this overestimation is 
predominantly a density effect (the gas expands due to repul- 
sive interactions and lowers the density, which in turn leads to 
a lower collision rate). At higher temperatures, the overesti- 
mation results from the density inhomogeneity of the trapped 
gas, which is not accounted for well by the scaling ansatz. We 
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FIG. 10. (Color online) Frequency and damping (top and bot- 
tom graphs respectively) of quadrupole oscillations in a 2D system 
(77 = 0) with N — 2200 particles. The thick colored lines are the nu- 
merical results for a 4th order basis set, including self-energy correc- 
tions. The thin black lines correspond to the analytic scaling ansatz 
analysis presented earlier (Sec.|V 2}. 



discuss this in more detail in the discussions section. 

Finally, we note that the same arguments apply to the quasi- 



2D case (Fig. 11 1 regarding the comparison between the scal- 
ing ansatz and the extended analysis. The additional fea- 
ture of the quasi-2D case, which was also discussed in de- 
tail in Sec. |V 2 is the reappearance of the collisionless limit 
at higher temperatures. A signature of this can be seen in 
Fig. [Tip by observing the non-monotonic behavior of the lo- 
cation of the peak of the damping rate. At higher tempera- 
tures, entrance to the crossover regime is delayed and the CL 
region expands. 



VII. DISCUSSIONS 

Most of the relevant discussions were already given in the 
main text. Here, we give a brief summary of the main re- 
sults and discussions, along with several complementary com- 
ments. 

In this paper, we started our analysis by investigating the 
equilibrium state of quasi-two-dimensional dipolar fermions 
in isotropic traps. In order to study the collective modes of the 
system, we solved the collisional Boltzmann-Vlasov equation 
for small perturbations of the trap potential with monopole 
and quadrupole symmetries. The self-energy corrections to 
quasiparticle dispersions and collisions were taken into ac- 
count via the self-consistent Hartree-Fock and Born approx- 



FIG. 11. (Color online) Frequency and damping of quadrupole oscil- 
lations for a quasi-2D system corresponding to r\ ~ 0.322 (refer to 
the caption of Fig.fTOlfor details) 



imations respectively. The validity of these approximations 



were assessed at the end of Sec. II C In particular, the usage 
of Born approximation restricts the validity domain to near- 



threshold scattering energies (see Eq. 25 1. We confined our at- 
tention to the regime where Tp <C ?dip, so that the scatterings 
remain in the near-threshold regime up to T ~ Tdip ^ Tp 
and make the thermal regime accessible to the scope of this 
work. 

We emphasize that once the conditions for the validity of 
CBV equation is met, this formalism is universally applica- 
ble to both CL and HD regimes, as well as the intermediate 
crossover regime. 

We carried out the analysis of collective modes in two 
stages: as a first approximation, we used bare quasiparticles 
and studied the response functions using the simple picture of 
scaling ansatz. This analysis implied that the monopole os- 
cillations occur at a fixed frequency of 2w > are undamped, 
and are independent of temperature and dipolar interaction 
strength. In case of quadrupole oscillations, however, we 
found a transition from the CL limit to the HD limit, indicated 
by oscillation frequencies of 2ui and respectively. We 

investigated the collisional relaxation rate of quadrupole os- 
cillations, v c , the single parameter that appears in the char- 
acteristic equation of quadrupole oscillations and yields the 
frequency and damping of the lowest lying quadrupole mode. 
This quantity was calculated for various temperatures and ver- 
tical trap frequencies and was shown to be expressible in terms 
of a universal function of T/Tp and 77. We found that in the 
2D limit (7/ = 0), v c is a monotonically increasing function 
of temperature and reaches to a plateau for large T/Tp. This 
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plateau persists up to T ~ Td; p beyond which the scattering 
enter the semiclassical regime and the scattering cross sec- 
tion starts to decrease as a function of temperature. The ex- 
istence of this plateau, which is a novel feature of dipolar in- 
teraction implies that (1) the character of trap excitations of a 
polarized 2D dipolar gas becomes weakly dependent on tem- 
perature in the regime Tp < T < Tdi P , and (2) collisional 
effects persists despite the fact that gas becomes very dilute. 
This behavior differentiates 2D dipolar fermionic gases from 
s-wave fermions where rarefaction of the gas at high temper- 
atures takes the system back to the collisionless regime for 
T > Tp. Also, the temperature window for collisional be- 
havior is universal for s-wave fermions and is not amenable 
to tuning, whereas for quasi-2D dipolar fermions, one can 
expand this window by (1) making the vertical confinement 
stronger to approach the 2D limit, and (2) either increase Td; p 
by using weaker dipoles or decrease Tp by using a weaker 
transverse trap. 

The existence of the plateau in v c is guaranteed as long 
as the scale separation Tp -C Tdip is met. By combining 
Eqs. ( f25] > and ( |68j ), one finds the condition for the plateau to 
lie in the collision dominated (hydrodynamic) regime: 

ATi < — < N*, (HD plateau) (82) 

ad 

The left and right hand sides of this inequality are equivalent 
to Tp <C Tdip an d N{ a d/ a o) 2 ^ 1 respectively, where 
the latter condition implies v c 1. The above inequality 
may be used as a simple experimental guideline to achieve 
hydrodynamics with dipolar fermions. 

In the second stage of calculations, we extended the anal- 
ysis by (1) including self-energy corrections and (2) satisfy- 
ing all of the moments of the CBV equation up to 8th order. 
Chiacchiera et al. ifTfl and Pantel et al. [37 1 have carried out 
similar extended moments analysis (up to 4th order moments) 
for the case s-wave fermions and have found the corrections 
to be significant and improve the matching between the theory 
and the experiment. 

This extended analysis allowed us to (1) investigate higher 
order modes for both monopole and quadrupole oscillations, 
and (2) evaluate the reliability of the simple scaling ansatz 
analysis. We found that despite the fact that satisfaction of 
higher order moments results in numerous new normal modes, 
the responses to the monopole and quadrupole perturbations 
(~ r 2 and x 2 — y 2 respectively) are predominantly governed 
by the lowest lying mode (in the quadrupole cases, we found 
contributions from a few overdamped modes as well). We re- 
mark that the frequency and damping of dominant mode, how- 
ever, is significantly modified by both self-energy corrections 
and inclusion of higher order moments. 

We argued that self-energy corrections is expected to re- 
sult in the damping of the lowest lying monopole mode, since 
the closure of moment equations that ensures the absence 
of damping relies sensitively on quadratic dispersions. We 
found that although this expectation is met, the damping re- 
mains very small (< 10~ 3 wo) even in the strongly interact- 
ing regime. The frequency of oscillations, however, is signif- 
icantly increased from its non-interacting value of 2wq. This 



correction was found to be most significant at lower tempera- 
tures where self-energy corrections are strongest. In order to 
rule out the possibility that the smallness of damping could 
be a result of the local effective mass approximation used in 
evaluating the collision integrals, we evaluated the collision 
integrals using the exact quasiparticle dispersion (altbeits, at 
the costs of a significantly increased computation time; see 
Appendix. |E|i for a few representative cases and found that the 
damping remains within the same order of magnitude. 

By investigating the velocity field of lowest lying monopole 
mode, we found that it retains its scaling character to an good 
approximation in all cases (i.e. v ~ r) with negligible temper- 
ature fluctuations. It is known from the hydrodynamic theory 
of non-ideal fluids that for a true isotropic and isothermal scal- 
ing flow, no dissipation results from shear viscosity or thermal 
conduction and the only source of dissipation is the bulk vis- 
cosity (for instance, see Ref. [38], §49). For such flows, one 
finds dS/dt = / d 2 rng ^"^(V • v) 2 where S is the total 
entropy and ( is the bulk viscosity. Note that the dissipation 
rate is small since it is second order in v. 

The extended analysis of the lowest lying quadrupole mode 
has the same qualitative behavior as the scaling ansatz analy- 
sis, albeit with significant quantitative corrections. At lower 
temperatures, self-energy corrections lowers the frequency of 
oscillations from the non-interacting value of 2wq proportion- 
ally to Xd- Also, we find that much stronger interaction is re- 
quired to reach the HD regime. This is simply a consequence 
of the rarefaction of the gas at low temperatures due to repul- 
sive interactions. Generally, we found that satisfying higher 
order moments results in a delayed entrance to the crossover 
(and the HD) regime. One may explain this finding by ob- 
serving that the scaling ansatz overestimates the collision rate 
since it describes the dynamics simply as a uniform time- 
dependent rescaling of the equilibrium distribution: in fact, a 
perturbation like x 2 — y 2 is most effective for the gas elements 
located furthest from the center of the trap while it does not 
affect the particles sitting at the center of the trap as strongly. 
Neglecting this fact clearly results in an overestimation of the 
collisional relaxation, specially since the gas is most dense at 
the center of the trap. 

In order to study the effects self-energy inclusion and 
higher moments separately, we have shown the quadrupole 
spectral function for three different temperatures and a fixed 
interaction strength using three different approximations in 
Fig. [12] 4th order basis set with and without self-energy, and 
1st basis set (scaling ansatz) with self-energy. At low tempera- 
tures (panel a), as one would expect, we find a significant cor- 
rection from the inclusion of self-energy. On the other hand, 
inclusion of higher order moments yields almost no correc- 
tion. At higher temperatures, this scenario is reversed. The 
reason is that self-energy corrections are essentially due to 
exchange interactions which diminish in the high temperature 
(classical) regime. On the other hand, the trapped gas assumes 
a Gaussian density profile at high temperatures (compared to a 
quadratic profile at T = 0) and naturally, higher moments are 
needed to describe the spatial inhomogeneity of the dynamics 
accurately. 

We also found that the most important corrections to the 
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FIG. 12. (Color online) A comparison of the quadrupole spectral 
functions obtained using three different approximations. The legend 
is shown on the plot (+£ and w/o E means with and without self- 
energy corrections). A^ = 0.45, 77 = and N = 2200 in all three 
plots, (a) T/T F = 0.1, (b) T/T F = 0.5, (c) T/T F = 1. 
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FIG. 13. Maximum relative deviations of the conserved quantities 
(N and E) for monopole oscillations for a sample configuration 
(T/T F = 0.1, A d = 0.5,r? = and N = 2200) 



scaling ansatz stems from the 4th order moments, beyond 
which the corrections become increasingly less significant. In 
practice, a second order basis set is sufficient to obtain the 
frequencies of the lowest lying collective modes within 0.1% 
range of the exact solution. Higher order modes naturally re- 
quire inclusion of higher order moments. 

As a consistency check for our numerical calculations, 
we investigated the satisfaction of the conservation laws 
(see Appendix [A|. The CBV equation conserves the particle 
number, momentum and energy, both in the differential form 
and the integral form. The quadrupole oscillations trivially 
satisfy these conservation laws due to the axial symmetry of 
the equilibrium state. This is not trivial for monopole modes 
since they have the same symmetry as the equilibrium state. 



Fig. 13 shows the maximum relative deviations of the particle 
number and energy in monopole oscillations as a function of 
moment satisfaction order for a sample case. We find that the 
particle number is conserved within a relative error of ~ 10~ 6 
even in a first order basis set (this is because one of the mo- 
ment equations is in fact the statement of conservation of 
particle number). On the other hand, we find that conservation 
of energy improves substantially upon expanding the ba- 
sis set. For the 4th order basis set, the relative error is ~ 10~ 5 . 



VIII. EXPERIMENTAL OUTLOOK 

The collective modes can be probed experimentally in var- 
ious ways. As described earlier, one common method is to 
perturb the trap potential with a short pulse and monitor the 
evolution of the cloud using in-situ or absorption imaging 
techniques (for example, see Ref. (9))- The relevant observ- 
ables are the radius and anisotropy of the cloud in case of 
isotropic and quadrupolar perturbations respectively. The fre- 
quency and damping of the collective modes are found by fit- 
ting the measured time evolution of the observable O exp (t) to 
a function of the form O m {t) = Aer^ sin(cjt) + Be" 700 ', 
where u> is the frequency of oscillations, and 7 and 70D are 
damping rate of the oscillatory and overdamped components. 
If required, the spectral function can be subsequently found by 
taking a Fourier transform of the measured impulse response 
signal cxp (t). 

Another method, which may yield more accurate results, is 
the measurement of spectral functions via modulation spec- 
troscopy. This is done by introducing a low-amplitude pe- 
riodic modulation to the trap potential at a fixed frequency 
il for a long duration r ^> c^ 1 , SI" 1 and measuring the 
absorbed energy. For a finite trap modulation pulse like 
SU ~ e~'*'/ T cos(fii) v(r), a simple linear response analysis 
yields 



AE : 



' abs ~ -T fl 3m[x v (r) + «/ r )] I 



(83) 



where Ai? a b s is the absorbed energy, v(r) is the shape of 



the trap perturbation (i.e. 



•j 2 and x 2 — y 2 for probing 



monopole and quadrupole modes respectively), and Xv(r) is 
the retarded correlator of v(r). Eq. (83i implies that the ab- 



sorbed energy in a modulation experiment yields a direct mea- 
surement of the spectral function. The absorbed energy can be 
measured in various ways. One method is to let the system re- 
thermalize after the modulation pulse, followed by mapping it 
to a non-interacting system by switching off the interactions 
adiabatically and finally measuring the rise in temperature of 
the non-interacting gas through a time of flight expansion ex- 
periment. The location of the peak in the spectral function 
and its width yield the frequency and damping of the collec- 
tive mode. 

We conclude this section by making predictions for the ex- 
periments with KRb. In the recent experiments with a quasi- 
2D configurations |4|, the vertical and transverse trap frequen- 
cies are uj z — (2ir) x 23 kHz and uj = (2tt) x 36 Hz re- 
spectively. The central layer has 2200 molecules, the tem- 
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perature is T = 500 nK and dipole moment is D = 0.158 
Debye, using which we find T/T F ps 4.36, r? w 0.322 and 
Arf s» 0.252. The dipolar temperature is T^p ~ 1.6 fiK and 
Tp/Tdip ~ 6.4 x 10~ 2 . Therefore, the near-threshold scatter- 
ing condition can be satisfied for quantum degenerate temper- 
atures. However, T/Ta p ~ 2.8 in the current experiments, 
implying that the system lies on the margin of semi-classical 
scatterings. Therefore, the predictions give here may not be 
very accurate. 

The mean-field corrections are small at this temperature 
and we can use the scaling ansatz analysis of Sec. V 2 
We find Q(f = 4.36, rj = 0.322) « 0.019 which yields 
v c w 0.04. Eq. (73 1 yields the frequency shift and damp- 
ing rate of the quadmpole mode as Slo/loq s» 1.25 x 10~ 4 
and 7 pa 0.01 uio ps 2.3 Hz. While the frequency shift is 
too small to be easily experimentally observable, the damping 
rate is sufficiently fast and can be easily observed. We remark 
that Q has a strong dependence on rj and rapidly decreases 
as f] is increased (see Fig. |5jl). Thus, one may increase the 
collisional relaxation rate significantly simply by making the 
vertical confinement stronger. In the hypothetical 2D limit, we 



,2D 



1.5 and consequently, 



.quad 



1.8 w « 415 Hz 



get v« 

and 7yi uad pa 0.3 wo ~ 71 Hz. The dramatic enhancement of 
the collisional effects in 2D compared to the quasi-2D con- 
figuration is remarkable. Also, we note that since T ~ Tdi P 
in the current experiments, the system is expected to remain 
in the plateau upon further cooling down to degeneracy at 
T~T F «115nK. 

Finally, we note that recent experimental progress with lan- 
thanide atoms such as 161 Dy |7| and 168 Er [3| which have 
large magnetic dipole moments (10 /xb and 1 jin respectively) 
is another possibility for studying collective many-body exci- 
tations of fermionic dipolar gases. 
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a consistency check for our numerical calculations. The con 
servation of mass (or equivalently, particle number) is: 

d f 

/ drn(p;r,i) = 0. 



dt 



(Al) 



The linearized equation, with the parametrization given by 



Eq. (39 1, yields: 



d 
dt 



y"drA $(p;r,t) -0. (A2) 

In the same parametrization, the conservation of momentum 
reads as: 



d_ 
d7 



J dr P A $(p;r,i) 



0. 



(A3) 



Finally, the energy density is given by £hf = p 2 /(2m) + 
mcjg/2 + E + [tt,]/2 in the Hartree-Fock approximation using 
which we get the following linearized form of the statement 
of conservation of energy: 

^ J dr (58 n + £ Ao$(p; r, t)) = °> ( A4 > 

where £ = Hq is the equilibrium energy density and 
58 = E + [5n}}/2 = S+[A $]/2. Using the properties of 
Hartree-Fock energy density functional, we get j dT58 tiq — 
(l/2)/drS+[A $]no = (l/2)/dTE+[no] A $. Using 
this identity, the two terms in Eq. ( |A4| l can be combined to 
yield: 



d 

dt 



drH A $(p;r,t) = 0. 



(A5) 



In case of quadmpole oscillations, these conservation laws 
as trivially satisfied due to difference between the symmetry 
of deviations and the equilibrium state. In the monopole case, 
while the conservation of momentum is still trivially satisfied, 
the mass and energy conservations may only be fulfilled if $ 
is a legitimate solution of the kinetic equation. 
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Appendix A: Conservation laws of the linearized collisional 
Boltzmann-Vlasov equation 



Appendix B: Asymptotic analysis of Q(T, 77 = 0) 

In the 2D limit (rj — 0), the asymptotic behavior of Q(T, rf) 
can be studied analytically. Setting rj = 0, the Erfcx functions 



appearing in the collision integral (see Eq. 67 1 evaluate to 1 



and the expression in the brackets in the second line simply 



becomes [\i — X2] 2 = sin £ sin v[l- 
will result in significant simplifications. 



sin 



. This 



1. Low temperature expansion 



The CBV equation admits three important differential con- 
servation laws for of mass density, mass current and energy, 
which can be established by multiplying the sides of CBV 
equation by 1, p and energy density 8 respectively, and inte- 
grating over p E0l . The collision integrals vanish identically 
in all three cases due to the existence of the same conservation 
laws in the level of 2-body scatterings. We state these conser- 
vation laws in their integral form here to utilize them later as 



In the low temperature regime, ji/T — > 00, we may use the 
following identity: 



lim {ft/ry 3 / P 5 dp 

1 



cosh(p — jtx/T) + cosh(fe 2/ o) 



1 

cosh(p — ji/T) + cosh(6ip) 
yiWW, (Bl) 
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in order to carry out the p integration. This identity can be 
established by observing that for large pp, the integrand will 
be exponentially unless p ~ ftp and 61,62 ~ {Pp)^ 1 . In 
the limit ftp — > 00, this implies that the right hand side must 
be proportional to <5(&i)#(&2). The proportionality constant 
can be found by integrating the left hand side over b\ and 62 
which yields the 47r 2 /3 prefactor. Identifying b\ and 62 as 
sin 2 ^ snv2v cos0 and sin 2 £ sin 2v cos(f>' respectively, we 
can carry out the £ and v integrations using the 5-functions 
and we finally get: 



Q{f ^ 0, 77 = 0) « C 



where C is given by: 



9 jo 



(B2) 



AM 2 



COS^ 



cos z 



and is equal to 19.176999 to six significant digits, 
be found analytically with little effort and we get: 



-8T 3 Li 3 (- 



(B3) 
)) can 



(B4) 



Using the asymptotic expansion of Li 3 (— x) for large x and 
the low temperature expansion of p mentioned after Eq. ( f33] >, 
we the following low temperature expansion: 



- Li 3 (-e^ f ) = 1/(6T 3 ) + tt 2 /(12T) + 0(T). 
Combining the last four equations, we finally get: 

Q(f -!■ 0,77 = 0)«;Cf 2 « 12.784666 T 2 , 



(B5) 



(B6) 



to leading order. This asymptotic limit is shown in Fig. [5^ as 
a blue dashed line and agrees with the numerical result. 



2. High temperature expansion 

The analysis of the classical limit (j3p — > 0) is simpler. 
First, we rewrite the hyperbolic functions in the denominator 
as cosh(p-mz) = e p /(2z) + (z/2)e^ p . Here, z = exp(p,/f) 
is the fugacity and goes to zero in the high temperature limit. 
Thus, cosh(p — laz) « e p /(2z) to leading order. In fact, 
the denominator of Eq. ( |67| is dominated by the first cosh 
term since the second ones are 0(e p ) -C e p / (2z), so that we 
can neglect them as well. With this simplification, all of the 
integrations become elementary and we get: 



Q(T^0,0) 



3tt)z 2 T 5 



(B7) 



Appendix C: Matrix elements of the evolution matrix in the 
monopole basis 



The variational linear response analysis of the CBV equa- 



tion, as described in Sec. IV requires the calculation of a large 
number of matrix elements. This task, however, can be some- 
what simplified as the angular integrations in the matrix ele- 
ments of M, 51 and Ho can be carried out analytically using the 
symmetries of the basis functions and the equilibrium state, 
reducing the problem to the evaluation of a two-dimensional 
integral over p and f . The latter computation can be done nu- 
merically accurately and efficiently. 

In this appendix, we provide readily computable formulas 
for the required matrix elements in the monopole basis. We 
define the shorthands R a = 2m a + k a , P a = 2n a + k a for 
given basis function <j> a . R a and P a count the powers of r and 
p in (j> a respectively. 



1. Matrix elements of M 



By definition, we have: 



1 Q /3 = / dT A (p, r) <f> a <j)p 

= J(2ir)fdfj^pdpA {p,f) f R * +R ? 

X pP«+ p n / (cosip) k " +kf < dip 
Jo 



E(k a + kf})(k a + k )\ 



2 



X p 



Aq(p, r) dr dp 



(CI) 



where E{n) — 1 for even n and E(n) = for odd n. For 
future reference, we define: 



h(n) = 



E(n) n\ 
2™ [(n/2)! 



(C2) 



The fugacity in the classical limit can be found from Eq. ( 33 1, an( j 
yielding z = 1/(2T 2 ) 
pansion — Li 3 (— z) = 



0(T 4 ). Using the asymptotic ex 
0(z 2 ), we finally get: 



Q(T ->oo,0)k -(3tt - 8) w 0.712389. 



(B8) 



This asymptotic limit is shown in Fig. [5^ as a red dashed line 
and is in agreement with the numerical result. 



C[A(p,f)] = J A(p 1 f)f m+1 p n+1 dfdp, (C3) 



using which we can write M a p — h(k a + hp) i, 



[Ao] 
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2. Matrix elements of Ho 

First, we evaluate the Poisson bracket {<pp, Ho}'- 

{0,3, Ho} = V f (pp • VpHo - Vp<^s • VrHo 
= 7p (P ■ V F )^ - 7r (r • Vp)0^ 

= 7p[2m /9 r/'(m fl -i,n fl ,fc jB +i) 

- 7 r [2n^j V'Cm^.n^-l.fe^ + l) 



(C4) 



where: 



7r 



r 2 r 



-2. 



7p = p- 2 p • VpH = 1 + p- 2 p • V P S . 
Plugging Eq. (C4 1 into the definition of (H ) Q( 3, we get: 



(C5) 



(Ho)a^ = y df A (j>a{^M 

= [2mp h(k a + kp + 1) + kp h(k a + hp — 1)] 
- [2ny3 fc(fc a + ^ + 1) + kp h(k a + kp- 1)] 



x I, 



(Pc+P/3-1) 



l/[7rA ]. 



(C6) 



3. Matrix elements of T 



By definition, 



E[A ^] = A c 



*(|p-p'U) A (p',f) Mp',f). 

(C7) 

It is easy to verify that a simultaneous rotation of r and p 
leaves E[Ao0^] invariant, so that Si?[Ao0^] may only depend 
on f, p and <p, the angle between r and p. Let cos ip = (p • 
p')/(pp r ) and cos <f> = (r-p)/(fp), so that r-p' = fp' cos(0+ 
ip). Expanding u(|p — p'|, rj) in a cosine series, 

00 

«(|p - P'l. V) = J2 u(n) (P,p'; V) cosM), (C8) 



n=0 



where: 



u {n) (P,P') = 



1 



7I"(^n,0 + 1) 



x J dip u ^ \/p 2 + p n — 2pp' cos ip, cosmp, (C9) 
and plugging into Eq. ( |C7| i, we get: 

S[A ^](p,f,</») = A d / ^A o ^,f)0 p "f R '> 



2tt 



00 /-27T j / 

V u(p,j?;r)) / — cos(m/0 cos(</> + ^) fe <*. (CIO) 
^0 ^ 2 ^ 



The angular integration can be evaluated using contour inte- 
gral techniques: 



(hp 
2^ 



cos(mp) cos(0 + ^) fc 



fc! 9(k-n)E(k + n) 



cos(n</>), (Cll) 



where 6(n) = 1 if n > and = otherwise. We denote 
the numerical prefactor in the brackets of the above equation 
by g(n, k). Plugging this into Eq. (CIO 1, we get: 



£[A <M(p,f» = Ad]T Qf ) fer)cos(n0), (C12) 



n=0 



where: 



dp' 



Q/T (ft f) = - 9 (n, fcp) J ^ A (p>, f) p 

xu^(p,ff;ri). (C13) 

The last integral can be easily evaluated numerically. Also, 
note that we only need v\ n > up to n = kp in order to evaluate 
f][Ao(pp] exactly. This is due to the fact that g(n, kp) vanishes 
for n > kp. Having evaluated E[Aq^], H a p can be evaluated 
readily by appealing to its definition: 

kp 

^ a p = Ad ^ ( [2m a g(n, k a + 1) + k a g{n, k a - 1)] 

n=0 

x [Qp n) A o7p ] - [2n a g(n, k a + l) (C14) 



+ k a g(n, k a - 1)] [Q ( ; ] A o7 r]) . 



4. Matrix elements of l c 



(C15) 



The evaluation of the matrix elements of the linearized col- 
lision integral operator is the most computationally expensive 
part of the calculation. In particular, the deviation of quasi- 
particle dispersion from the bare quadratic dispersion makes 
the calculations even more challenging. To our knowledge, 
all of the previous works along this line have evaluated the 
collision integrals for bare particles. This approximation is 
justified when one is dealing with the the Boltzmann equation 
where one neglects mean-field corrections altogether. How- 
ever, since we have included mean-field effects on the dynam- 
ics, we must also use the dressed quasiparticles dispersion in 
order to satisfy conservation of energy. In order to do this 
in a numerically tractable way, we have found that the quasi- 
particle dispersions can be approximated well using a local 
effective mass approximation (within an error of less than 2 
percents). To this end, we approximate the dressed quasipar- 
ticle energies as: 



H {p,r) sa e (r) 



P 



2m* (r) 



(C16) 
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where: 



e (r) = S o (f;0), 



m*(r) 



l + d 2 p Z (f;p) 



-i -l 



p=0 



(C17) 



As we will see shortly, this approximation allows us to put 
the collision integral into a simple form suitable for numerical 
treatments. As a first step, we go to the center of mass frame 
of the colliding particles and define: 



P= 2 



p' 

P' = y + q' 



Pl = 2 q ' 
P' 

pi = Y -q' 



(C18) 



using which we get: 



d 2 r 



d 2 p d 2 p! d 2 p' d 2 p[ 

(2tt) 2 (2tt) 2 (2tt) 2 (2^) 2 



(2ir)6(AE) (2tt) 2 (5(AP) 



m*(f) PdP qdqd^d^ 

— ^ — rdrdip— — — — — , (C19) 

2 y 2tt 2tt 2tt 2tt ' 

where 0, 0' and ip are defined as cos <j> = q-P/(qP), cos </>' = 
q' • P/(q'P), and cos ip = f • P/(fP). Note that P = P' and 



<7 = g' in the remainder of the integrand due to conservation of 
momentum and energy respectively. The scattering amplitude 

M = XdH\p-p'\,v)-u(\p-p[\,v)] -> A d [w(2g| sin[(</>- 
0')/2]|,jy) - u(2q\ cos[(<j> - ^')/2]|,7?)]. The product of the 



equilibrium distribution functions, n no,i(l 
can be conveniently written as: 



n a n ,i(l - rioX 1 - «o,i) 

1 1 



«o)(l 



4 cosh E + cosh 7 cosh E + cosh 7' 



where E = /3(P 2 /4 + q 2 ) /[2m* (f)} + /3f 2 /2 - /3fl, 7 = 
/3P<7 cos <p/[2m* (f)], 7' = f3Pqcos ^>'/[2m*(f)]. The an- 
gle V is only present in S [</>«] S[</>J. Therefore, the in- 
tegration over if) is immediate and elementary. We eval- 
uate the required integrals using Mathematica and define 
S a p(f,P,q,^>,^)') = J dipS[(j) a ]S[(j)p]. The integral can 
be put in a more useful form using the change of variables 
P = (8/?//?) 1 / 2 sin£ cos^, q = (2p/^) 1 / 2 sin£ sin;/ and 
f = (2p//3) x l 2 cos£, where p e [0, 00), 1/ e [0,tt/2] and 
£ G [0, 7r/2]. The final expression is: 



'a/3 = 



(2iV)H 2 



/*00 /'2'7r i 1 pIlT 1 if p ^ p ^ 

/ P 2d P / ir- IT- d£sin 3 £cos£ / dv sin2f 
Jo Jo 2?r Jo 2?r Jo Jo 



8(2 7 r) 2 /3 Ar °+^+ 3 
x S Q( 3(i/2pcos^, -yS/Jsin^ cos z/, -\/2psin£ sin;/, 0, </>') m*(f) 



0u y2^/2p/P sin£ sinz/|sin[(0-^)/2]|,r?J - \Jf3u \2^2p//3 sin£ sin 1/ | cos[(<f> - <j>')/2]\, t] 
x { [cosh (psin 2 £/ra*(f) + pcos 2 £ + f3eo(f) — Pp) + cosh (p sin 2 £ sin 2v cos (f>/m*(f))] x (</> o </>')} 



(C20) 



where TV, 



a{b) 



^a(6) 



la(I>) 



fc a(f) ) and r 



in m*(f) and e(f). We evaluate the above 5-dimensional in- 
tegral for all pairwise combination of basis functions using a 
numerical Monte-Carlo integration with 5 x 10 8 points which 
we found to yield a relative statistical error of less than 1CT 3 
in all cases. 



v / 2pA3cos£ (jia.iA,) = (0,2). 



1. Matrix elements of M 



Appendix D: Matrix elements of the evolution matrix in the 
quadrupole basis 

In this appendix, we provide readily computable expres- 
sions for various matrix elements in the quadrupole basis by 
carrying out the angular integrations analytically. For a given 
quadrupole basis function £i</> a , we define a pair of num- 
bers (/Xj,fj) as the number of powers of r and p in £3 re- 
spectively, i.e. (01,^1) = (2,0), (02,^2) = (1,1). and 



The angular integrations in M can be easily carried out us- 
ing the variables cos^ = f • x and costp = f ■ p/(rp)- In 
this variables, we get £ 8 = f IXi p Vi cos(2(f> + Vjip). The angular 
integration are elementary and we find: 



dTAo&&. 



Vj I , k a 



(Pc+Pfi+Vi+Vj) [ 



(Dl) 
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2. Matrix elements of H 

As a first step, we evaluate the Poisson bracket 
{£j(t>f>,Uo} = ^j{<l>p,'Ho} + 4>p{^,'Ho}- The expression for 
{0/3, Ho} is known from the previous appendix (Eq. C4i. We 



find {£j,'Ho} = Xjk{p, r) (sum over fc is implied), where: 



Xjk — 



I 2 7p \ 
V -2 7r / 



(D2) 



Thus, we get: 



(Ho)^ = J drA o ^0 Q {^^,Ho} 

df A <M</> / 3,'HoKi£.? 



drA , 



(D3) 

The angular integrations in (Ho)^ j can be most easily eval- 
uated using the variables defined earlier, cos <fr — f • x and 
cos ip = f • p/(rp). The final result is: 

( h o)q/j,i = ^ 2mpg{\v i -v j \,k a + kp + l) 
+ kp g{\v, -Vj\,k a +kf)- 1) I^+Pe+vi+^i) t A o 7*>] 
2npg(\i/i - fjljfca, + fc/j + 1) 



+ fc^ g{\vi -Vj |, fc a + kp - 1) 



7 -(fl Q +fl. /3 +M 1 +^ + l)r A ^ 1 
1 (P a +Pp+u i + Vj -l) i A 7rJ- 



The angular integrations in (Hq)qo 2 are similar 10 those in 
(M)^. The result is: 

( H o)^, 2 = g 5(1^ ~ v k\,K + k ) 



where we have expressed u(|p — p'\,rj) in a cosine se- 
ries. The angular integration (the second integral) can 
be carried out using the contour technique and yields 

g(vj, n, kp) cos(2</>) cos(ra/>) — h(fj, n, hp) sin(2</>) sin(m/>)> 
where: 

t7(0,rc, fc) = g(n,k), 

g(l,n,k) = g(n,k + l), 

g{2,n,k) = 2g(n,k + 2)-g(n,k), 

h(v, n, k) = g(i>, n, k) — g(y + n, k). (D6) 

Plugging this back into Eq. {D5}, we get: 

S[A O/j] = Qp n j(p,f)\g(vj,n,kp)cos(2<j>) 



71=0 



x cos(n-i/>) — h(fj, n, kp) sin(2</>) sm(nip) 



, (D7) 



where 



xu(p,p';ri)A (p',f). (D8) 



The last integral can be evaluated easily numerically. The final 
result can be expressed easily in terms of the last two expres- 
sions: 



n=0 



1 n (u k ,n,k a ) r (R a +Vk) r A r>( n ) Y" 1 



2 W I-',.".A .1 "" V V 



- ( 2m Q g^.".^^) + k a G^'^-V) 



v r(iJa+/i«+l) [A nWn/ 



(n). 



, (D9) 



where we have defined the shorthand notation (j , ) I/1 ' ni, J0 
g(isi,n 1 ,k 1 )g(v 2 ,n2,k 2 ) + h(iy 1 ,n 1 ,k 1 )h(iy 2 ,n 2 ,k 2 ). 



3. Matrix elements of T 

Like the monopole case, the first step is evaluating 

S[Ao6^»]: 



A 



2tt 



x / ^ cos fc "(V> + ip') cos[2^ + i/,-(i/j + ip')] cos(ntp'), 
' 2tt 



(D5) 



4. Matrix elements of l c 

The matrix elements of the collision integral in the 
quadrupole basis is identical in form to those of the monopole 
basis (Eq. C2Q\ . The only differences are (1): S a p must be 
replaced with: 

/c\f) 
~d^Sfe<ySfo<M, (D10) 

where we introduced an extra angle cos# = e x ■ P / P. (2): 
The prefactor ^ a +N ff +3 pN a +N ff +5 m the denominator 
due to the extra powers of and £j introduce. The defi- 

nition of N a (p-j is the same as before. 
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Appendix E: Calculation of the collision integrals with exact 
Hartree-Fock quasiparticle dispersions 




FIG. 14. (Color online) The damping rate of the monopole oscilla- 
tions in 2D and with N = 2200 particles, (a) T/T F = 0.5, (b) 
T/Tf = 1.0 and (c) T/T F = 1.5. The thick colored lines are the 
previously discussed result obtained using the local effective mass 
approximation of the the dispersions in the collision integral. The 
dashed lines denote approximate solutions obtained by relaxing the 
conservation of energy (from top to bottom, a — 0.05, 0.02, 0.01 
and 0.005). The thick black line is extrapolation to a — (the exact 
solution). 



In Sec. |C 4| we simplified the expression for the collision 
integral matrix elements using the local effective mass ap- 
proximation (LEMA) for the quasiparticle dispersions. Al- 
though we found this scheme to be a decent approximation 
in the weakly interacting regime (the approximate dispersions 
lie within a few percents of the exact Hartree-Fock disper- 
sions), one may argue that an exact treatment is necessary for 
stronger interactions. This objection is more serious when one 
is looking at the effects that crucially depend on self-energy 
corrections, such as the damping of the monopole mode. In 
this section, we address this issue and present numerical jus- 
tification for the reliability of LEMA. 



The major simplification resulting from LEMA is the pos- 
sibility of an analytic treatment of the 5-function in the col- 
lision integral associated to the conservation of energy (see 
Eq. C19 1. In that case, one simply finds q = q' , where q and q' 



are the magnitude of the momenta of the initial and final scat- 
tering pairs in the center of mass frame. Without a spatially 
local quadratic dispersion, this result does not hold anymore 
and in general, there is no easy way of treating the 5-function 
analytically. Here, we adopt a simple numerical approach to 
overcome this difficulty. Using a limiting process to to define 
the delta functions, 



S(AE) 



lim 

(T-S.0 



1 



-AE 2 /(2a 2 ) 



2na 



(El) 



we may replace the <5-function with Gaussians and calculate 
the collision integrals for various values of a. The a — > 
limit may be found by extrapolating the obtained results. This 
approach is considerably more computationally demanding 
than LEMA, however, it yields an accurate calculation of 
the collision integral matrix elements. The integrals are six 
dimensional in this case (over the variables f, P, q, q', <f> and 
(/)') since q and q' may have different values now. 

We implemented the above method for the case of 
monopole oscillations for a 2nd order basis set (satisfying all 
of the 2nd and 4th order moments of the CBV equation). The 
extrapolation is carried out using a polynomial fit. Fig. [T4| 
shows the damping of monopole oscillations obtained using 
several choices of a, the extrapolated result, and the result 
obtained from the effective mass approximation (shown ear- 
lier in Fig. |6j. The matching between the effective mass ap- 
proximation and the exact result is excellent up to A <j ~ L 
The LEMA result, however, deviates from the exact result for 
Xd > 1. In any case, we find 7^0^' < 10 _3 wo and our con- 
clusion about the smallness of the damping of the lowest lying 
monopole mode is still valid. Finally, we note that the im- 
provement in the prediction for the frequency of oscillations 
is much smaller (a relative correction of about 10~ 6 ) even in 
the strongly interacting cases. This is due to the fact that the 
frequency shift arises essentially from the self-energy correc- 
tions on the dynamical side of the CBV equation, which is 
already treated exactly. 
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